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Diffusion of Cobalt in KCI Single Crystals 


DOMINIC DAMIEN & K V REDDY* 
Department of Physics, Indian Institute of Technology, Madras 600036 


Received 19 June 1982; revised received 2 September 1982 


Results of the studies on the diffusion of cobalt in KCI are reported. Annealing was restricted to the temperature range 700- 
760°C, due to the poor solubility of cobalt. While profiles of crystals annealed for short durations exhibit only one region, long- 
time annealed samples show up a region of slow diffusion near the surface, in addition to the region of normal diffusion. This 
anomalous character is explained as due to the poor solubility of CoCl3~ ions formed in alkali halide lattices. The region of 
normal diffusion is treated as due to the migration of this complex ion by vacancy mechanism. Values of the activation energy 
for CoClj~ ion diffusion and the enthalpy of migration, obtained from the present study, are 1.71 eV and 1.12 eV respectively. 


1 Introduction 

Alkali halide crystals have been studied extensively 
to understand the existence of various types of defects 
in monocrystals. Often alkali halide crystals are doped 
with divalent cation or anion impurities to determine 
the nature of point defects present. It has been 
established''? that divalent alkaline earth impurities 
such as Ca?*, Ba?* and Sr?* enter alkali halide 
lattices substituting normal cations. When a divalent 
cation replaces a monovalent host cation, a cation 
vacancy is created in order to satisfy charge neutrality 
of the crystal. Divalent anion impurities such as SO2Z- 
and CO2~ occupy anion sites leading to creation of 
charge compensating anion vacancies*’*. 

It was speculated that unlike alkaline earth 
impurities, cobalt enters NaCl and KC lattices as 
interstitial CO?* in view of its small ionic size (ionic 
radius, 0.72 A) as indicated from earlier optical 

: absorption studies*’®. However, dielectric loss 
"measurements on cobalt-doped KCI [Ref. 7] and NaCl 
[Ref 8] have been interpreted assuming Co*” as 
a _ substitutional impurity. Allnatt and Pantelis’ have 
studied diffusion of cobalt in NaCl crystals and found 

that the diffusion profiles are qualitatively different 

from those of other divalent impurities. These profiles 
y a region of anomalously low diffusion near the 
e followed by a second region of normal 
yn. They have not given any explanation for this 
yrdinary behaviour of cobalt. Nor could they 


Y. 


whether Co** occupies substitutional site or 


“the workers have ignored the fact that 
form a complex ion CoCli~ in KCI and 


a ae 


similar to SO2~ and CO, and can enter halide 
lattices substitutionally replacing a host anion. Our 
ionic conductivity measurements’! on cobalt-doped 
KCI have confirmed the creation of charge 
compensating anion vacancies. In this paper, we are 
reporting the results of our studies on the diffusion of 
cobalt in KCI. The analysis has been carried out taking 
into account the formation of the complex ion, 
CoClz~. Due to poor solubility of cobalt in KCI [Ref. 
11], we had to confine the diffusion studies to the range 
700-760°C. 


2 Experimental Details 

Single crystals of KCl were grown from melt by 
Kyropoulos technique starting from AnalaR grade 
material. The accidental divalent impurity con- 
centration was estimated to be in the range 0.5-1.5 
x 10~°, from electrical conductivity measurements. 
This concentration of divalent impurities does not 
have any significant effect on our results, since all the 
diffusion studies were done above 700°C. Crystals were 
cleaved into blocks of about 1 cm x 1 cm x 0.3 cm size 
and one of the faces of each of these samples was 
microtomed. A few drops of CoCl, solution, 
containing radioactive®° Co as tracer, were transferred 
onto the microtomed face. The solution was allowed to 
evaporate to dryness, while a thin layer of radioactive 
CoCl, of about 2m thickness was formed on the 
crystal surface. Two such samples were put together in 
a graphite box, wrapped in a platinum foil with their 
active faces touching each other. 

Annealing was done in argon atmosphere by 
inserting the graphite box into a quartz tube. Thermal 
gradient across the sample was minimized by 


introducing the quartz tube into a stainless steel cavity 
as shown in Fig. 1. The annealing temperature was 


maintained within +0.2°C employing an electronic 
PID (Proportional, Integral and Derivative) tempera- 
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Fig. |—Set-up for diffusion annealing in argon atmosphere (Figure 

not to scale) [GB, graphite box; QT, quartz tube filled with argon: 

SB, stainless steel block to reduce thermal gradient; TC, 
thermocouple; FC, heater coil; TI, thermal insulation ] 


ture controller. (The controller has been designed and 
constructed in our laboratory using operational 
amplifiers and silicon-controlled rectifier.) A calibrated 
chromelalumel thermocouple was used to make 
temperature measurements and the temperature was 
recorded continuously. After annealing for sufficient 
time, the crystals were quenched to room temperature. 
The four sides of the active face were cleaved off to 
avoid the effect of surface diffusion. 

The crystals were then microtomed (Spencer 
Microtome, Model SI-750A) and the material removed 
was collected in flat dishes. To reduce self-absorption 
of beta particles, the collected material was dissolved 
by adding a few drops of water. On drying up, the 
material was left as a thin layer. Activity in each dish 
was determined using a GM counter. Diffusion profiles 
were plotted in the usual manner—logarithm of 
activity per micron as a function of square of the 
distance from the surface. 


3 Results 

The profiles obtained for crystals, which were 
annealed for a longer time, were all similar to what 
Allnatt and Pantelis® have Observed for cobalt 
diffuison in NaCl. As can be seen from Fig. 2, the 
Profiles have two distinct regions. The region near the 
surface has a very high slope indicating very slow 
diffusion. The value of diffusion coefficient for that 
region turns out to be of the order of 10-12 cm?/s. For 
the second region, the value of diffusion coefficient is of 
the order of 10~'° em?/s. The experimental points lie 
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Fig. 2—Diffusion profile of cobalt in KCI after annealing for 165 hr 
at 752 C [Fitted curve is shown by the solid line.] 
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Fig. 3— Diffusion profile of cobalt in KCI after annealing for 76 hr at 
721°C. 


scattered for this region mainly due to rather low count 
rate. 

Crystals which were annealed for shorter durations 
gave profiles that appear quite different. If we exclude 
the highly active first slice, the Profiles have only one 
region (Fig. 3). The value of diffusion coefficient is of the 


order of 10°'° om?/s and it indicates that the first 
region is missing. : 


DAMIEN & REDDY: DIFFUSION OF COBALT IN KCI CRYSTALS 


4 Analysis 

As mentioned in the introduction, cobalt has a 
tendency to form the complex ion CoCl2~. The 
chemical reaction equation can be written as 


CoCl, + 2KCI> K,CoCl, 


The solubility of CoCl2~ ions is very poor in KCl 
lattice'', which is a common property exhibited by 
divalent anion impurities*. Excess of cobalt may be 
present in precipitate form of composition K ,CoCl,. 
The dissolved CoCl{~ ions occupy anion sites creating 
charge compensating anion vacancies. This ion can 
diffuse into KCI just like SO3” or CO3™ ion, by 
vacancy mechanism. And the second region of the 
diffusion profiles (Fig. 2) corresponds to the diffusion of 
this divalent anion. Near the surface, cobalt may be 
present as CoCl, and K ,CoCl,. CoCl, can diffuse and 
then react with KCI (solid-solid diffusion reaction) to 
form K,CoCl,. It is also probable that the reaction 
may be taking place at the surface itself and then the 
molecule, K,CoCl, (in precipitate form) diffusing as 
such via dislocations. However, the first region can be 
treated as the reaction profile. 

We have tried to fit the experimental profiles with 
the expression : 


A(x) = A, exp(— x?/4D,t) + A,exp(— x?/4D,1) 


for activity at a depth x from the surface. D, and D, are 
the diffusion coefficients for the first and second regions 
and f¢ is the time of annealing. The fitted curves are 
shown as solid lines in Figs 2 and 3. Fitting is quite 
_ good for the first region, while the points lie scattered 
for the second region (Fig. 2). This is as expected, since 
the count rate is rather low for the second region. Short 
ses anneal-data also could be fitted with the same 
ssion, though there are no experimental points in 
irst region other than the first point (Fig. 3). This 
1ins why short time-anneal profiles do not exhibit 
zion of slow diffusion. 

nius plots of diffusion coefficients for both the 
ire shown in Fig. 4, for the temperature range 


wo 


ra for divalent ions (10-2-10~* 
Bares panes to A vel 
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Fig. 4— Arrhenius plot of the diffusion coefficients of cobalt in KCI 
[Curve: |. for the first region of slow diffusion: and 2. for the second 
region of normal diffusion] 


eV)’. It is even higher than the value reported for SO,” 
diffusion (1.23 eV) in KCI [Ref. 3]. 

As proposed already, the second region of diffusion 
profiles has been treated as due to the migration of 
CoCl{ ion by vacancy mechanism. Diffusion of 
divalent impurities is by the exchange of impurity and 
the associated vacancy and hence depends on the 
degree of association (p) of the impurity with vacancies. 
It also depends on the jump frequency of a tracer (@) 
into an associated vacancy and the jump frequency ofa 
lattice ion (@,), which is nearest neighbour to both a 
tracer and the associated vacancy, into the vacancy. In 
the intrinsic range of temperature, D is independent of 
pam impurity concentration and can be written as 

D =(a’/3)o3 fp 
where a is the anion-cation separation and f is the — 
correlation factor. Since @,>,, the correlation 
factor tends to unity. 

Temperature dependence of p can be written as 


es Xo K3/(1 + X90 K3) 

where the concentration (Xo) of free anion vacancies is 
given by, 

Xo = exp (Ss/2k).exp(—hg/2kT) 


and 
K, = 12exp(s,/k).exp(—h,/kT) 
eens dependence of impurity jump frequency 


0 = exp(- se 
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Taking the values for entropies and enthalpies from 
our conductivity studies'! (ss/k = 8.00, hs = 2.49 eV, 
s,/k=0 and h,= —90.54 eV), we have fitted the 
diffusion coefficients as a function of temperature by 
the method of least squares, treating @o and Am as the 
parameters to be determined. The best fit is given by, 


w, = 7.66 x 10 exp(— 1.12/kT) 


The value of enthalpy of migration turns Out 2 be 
higher than that for divalent cations (0.5-0.9 eV)”. 


5 Discussion 

Our results are similar to the observations of Allnatt 
and Pantelis? on the diffusion of cobalt in NaCl. Iida 
and Tomono!? did not observe the first region 
probably due to short annealing time as pointed out 
earlier. In short, the first region corresponds to the 
formation of CoClq ions and the second region 
corresponds to the diffusion of these ions. 

The present study favours the proposition that 
cobalt forms a tetrahedral complex ion (CoCI2~) in 
‘KCI crystals. The poor solubility of this ion is 
responsible for the two distinct slopes of the diffusion 
profiles. Had the solubility been higher, the first region 
would have been submerged in the second region. In 
fact, a region of slow diffusion is not observed for 
calcium and strontium, which are highly soluble in 
halide lattices. . 

Beniere er al.° investigated the diffusion of SO27 in 
KCI and NaCl and obtained vaiues for diffusion 
coefficients higher than those for self-diffusion 
coefficients. The higher values have been accounted for 
as due to the effect of the large size of the ion on the pre- 

’ exponential factor, which is a general behaviour 
reported for diffusion of monovalent ions!®, We also 
have observed the pre-exponential factor higher as 
expected. But the high value for the enthalpy of 
a (1.12 eV) makes the diffusion coefficients 


"1 
vale - 7 


“a The Present study, by itself, does not prove that 
E = alt forms complex CoCl2~ ions and that these 


e * / 


er wre hah & 


radicals occupy anion sites. But, when we combine the 
results with those of the electrical conductivity studies 
of cobalt-doped KCI [Ref. 11], we find some 
convincing evidence in this regard. The poor solubility 
cannot be explained, if one assumes that Co?* 
occupies cation site or interstitial position. Further 
evidence can be had from the measurement of self- 
diffusion coefficients of cobalt-doped crystals in the 
extrinsic region. If CoClq” ions occupy anion sites, 
charge compensating anion vacancies are created. This 
will enhance the diffusion coefficient of Cl” and 
suppress the diffusion coefficient of K* in KCI crystals. 
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Magnetic Properties of Divalent Nickel Ions in Lead Silicate Glasses 


A L HUSSEIN, F A MOUSTAFFA*, Z ABD EL-HADI** & A G MOSTAFAt 
Glass Research Laboratory, National Research Centre. Dokki. Cairo. Egypt 


Received 29th May 1982 


Magnetic measurements of some lead silicate glasses containing nickel oxide were carried out to study the equilibrium 
between the octahedral and the tetrahedral coordinations of the divalent nickel ions. The effects of the composition of the base 
glass and of the concentration of the nickel oxide were studied. The results indicated that these glasses are paramagnetic and 
their paramagnetic susceptibilities in the octahedral coordination are lower than those in tetrahedral coordination. The 
equilibrium between these two states of the nickel ion coordination depends on several factors such as the polarizability (of the 


oxygen ligands surrounding the nickel ion), 


mobility (of the alkali ions), field strength and the ability of some divalent, 


trivalent and tetravalent metal oxides as network forming units (MO, or MQ,). 


1 Introduction 
The magnetic susceptibility measurements of 
coloured glasses are of particular interest to follow the 
different states of coordination in which the elements of 
transition-metal group can exist in glasses. The 
influence of the chemical combination on the magnetic 
properties of compounds and on the compounds of the 
transition elements was summarized by Selwood’. 
Studies on paramagnetism of the transition elements in 
glasses?~° have been reported and Breit and Juza’ 
measured the magnetic properties of Ni in glass. 
Moore and Winkelmann* assumed that nickel exists in 
glass as divalent Ni ion in association with oxygen 
atoms. Juza and Schultz® studied the magnetic 
susceptibilities of Ni in different types of glasses, and 
attributed the colour changes to the change in the 
coordination number of Ni. The magnetic properties of 
Ni ions in glasses of various compositions have been 
explained on the basis of the ligand field theory”. 
= Studies comparing the measured susceptibility values 
with the theoretically calculated ones for the different 
- jons in glasses were reported by several wor- 
a s2~5:719 The magnetic properties of the divalent 
‘Ni ions in different types of glasses such as alkali 
yrate, lead borate and lead silicate have already been 
, by one of the present authors'?. In this paper, 
netic properties of divalent Ni ions in lead 
glasses are studied and the effects of 
of the base glass and concentration of Ni 


scussed. 


) tele 


si tal Details and Calculations 
pure grade materials were used for 
1l samples of glasses. The samples were 


University College for Girls, Ain-Shams 
a ae 


gYp' 
7 


finally pulverized and then melted in Pt-2%Rh 
crucibles in an electrically heated furnace at a 
temperature of 1100-1200°C, depending on the glass 
composition, for 4 hr. The molten samples were then 
poured as rods (1 cm x | cm x 10cm) and after setting, 
they were annealed. The glass samples were ground 
and polished in the usual way with minimum amount 
of water. The magnetic susceptibility of these glass 
samples was determined following the Gouy method 
using an electromagnet of approximately 14700 G field 
strength. The volume susceptibility (K ) for each sample 
was determined using the equation'?: 


K = K' + 2Amg/aH? sf 


where K’ is the volume susceptibility of air 
(approximately zero); Am the weight difference due to 


the magnetic pull; g the acceleration due to gravity; a — 
the cross-sectional area of the glass sample and H the > 


magnetic field strength. Then the value of the mass 
susceptibility (~) for each sample was calculated using 
the equation : 

y =K/p .--(2) 
where p is the density. Similar measurements were also 


carried out on the base glass (free from Ni) for each 


composition, and the difference gives @ values for Ni 
ions. 


3 Results 

The composition of the first set of samples ranged 
from PbO 70%, SiO, 30% to PbO 90%, SiO, 10%, 
each glass sample containing 0.25 g nickel oxide/100 g 
glass. Our results indicated that all these glasses are 
paramagnetic and @ increases with the gradual 
increase of PbO content (Fig. 1). 

The samples containing different concentrations of 
NiO indicated that all these glasses are also 
paramagnetic. The values of ¢ increased as the NiO 
content was gradually increased as shown in Fig. 2. 
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Fig. 2—Plot of @ versus NiO content in the glass samples studied 


Mass Suscepx1o 


R,0 % 


Fig. 3—Effect of replacement of PbO by alkali oxide on the value of 
@ in the glass samples studied 


When lithia, soda or potash was introduced 
replacing 1, 5 or 10 parts of PbO, the values of 
increases, and the relative magnitude of the effect in the 
case of various oxides was in the order: Li,0 > Na,O 
> K,0 (Fig. 3). : 

When magnesia, zinc oxide, lime, strontium oxide or 
barium oxide was introduced replacing 1, 5 or 10 parts 
of PbO, the value of @ increases, and the order of the 
relative effect was: MgO > ZnO > CaO > SrO > BaO 
(Fig. 4). 

When alumina was introduced replacing | or S parts 
of silica, and when titania or zirconia was also 
introduced replacing 1 or 1/10 Parts of silica 
respectively. the @ values changed as shown in Table I 

Generally, the colour of all the Studied samples was 
reddish brown: as the nickel Oxide content was 
increased, the colour changed gradually from faint 
reddish brown to deep reddish brown. 


"0 
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A 
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Fig. 4—Same as in Fig. 3 but for Ca, Sr. Ba, Mg and Zn oxide 
replacements 


Table 1!—Compositions and Values of Magnetic 
Susceptibilities (gy) of Lead Silicate Glasses 
PbO SiO, R,O, or RO, @ x 10° 
80 19 1 (AI,O3) 0.1778 
80 15 5 (Al,O,) 0.2317 
80 19 1 (T10) 0.1522 
80 19.9 0.1 (ZrO;) 0.1482 


4 Discussion 

The paramagnetism introduced into a glass by the 
incorporation of the transition-metal ions arises from 
the unpaired electrons of the unfilled 3d-subshell of 
these metals. The d-shell contains five orbitals, each of 
which may be occupied by one or two electrons. The 
distribution of electrons among the five orbitals is 
governed by Hund’s rule for maximum multiplicity to 
produce the maximum possible spin and angular 
momentum which is achieved by having the minimum 
number of paired electrons!2. 

Nickel was found to be paramagnetic in either its 
atomic or ionic state and it has paramagnetic 
properties when it occupies the octahedral or the 
tetrahedral coordination state*. On this basis, the 
experimental results obtained for our glass samples can 
be discussed as follows, 

In the glasses containing low proportions of PbO, 
the Pb ions are enclosed within the interstices formed 
by SiO, tetrahedra. Thus, it can be expected that Ni 
will be found in such glasses in the octahedral 
coordination state. As the PhO content was increased, 


4 considerable proportion of Pb ions will act as double 
bridges between two adj 


t SiO tetrahedra . 
groups, as = Si~O—Pb—O—§} a may be paar 


‘4 2 ae 
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besides the formation of the PbO, groups. Since the Pb 
ions are easily polarizable, the structure becomes very 
Open, easing the movement of the atomic aggregates. 
These considerations may account for the tendency of 
the Ni ions to be present in the tetrahedral 
coordination state in the glasses containing high 
proportions of PbO. This interpretation was found to 
be in agreement with the experimental results obtained 
(Fig. 1). 

The increase in the values of g as the NiO content 
was increased should be attributed to the para- 
magnetic character of the divalent nickel ion Ni? * itself 
due to the presence of the unpaired electrons. 

When soda, lithia or potash are introduced in the 
glasses, these cations can occupy the holes in the glass 
structure which cause a change in the forces exerted on 
the oxygen ions of the SiO, groups. Not all the Si** 
ions will be directly linked to each other, as in silicate 
glasses, but will become linked together through PbO, 
or bridges of Pb** ions. The flow of lead silicate glasses 
would be expected to be controlled by the silicon- 
oxygen bond. The presence of alkali ions (in such glass) 
would weaken these bonds: (i) because the silicon 
atoms are no longer in the field of tetrahedral 
symmetry and (ii) the alkali ion has a polarizing effect 
on the oxygen around Si** ions. 

| Thus, on replacing Pb ions by one of the alkali ions 
-__ (R*), the polarizability of the oxygen ions will decrease. 
Accordingly, the proportions of the Ni ions in the 
-_ gctahedral coordination will decrease. The relative 
effect was in the order: Li,O > Na,O > K,0O and it 
depends on the polarizing power of these ions which 
es exhibit on the neighbouring oxygen ligands and 
. | the compactness of the alkali cations in the 
her 1 magnesia, zinc oxide, lime, strontium oxide or 

m oxide are introduced in these glasses, these 
nl t cations may form bridges in the glass structure 
' ey be enclosed within the interstices in the 
work. Magnesia and zinc oxides can partially 
lass forming groups like MgO, or ZnOg. 
placement of PbO by one of these divalent 
25 would cause a decrease in the BO,/BO, 
lingly the proportion of the octahedral 
ckel ions will increase with the gradual 


° 


increase of one of the divalent metal oxides replacing 
PbO. The effect of these different divalent metal oxides 
would be expected to depend on their coordination 
State, their concentration and the addition of other 
structural groups. The relative effect of these oxides 
was as follows: MgO(highest) > ZnO > CaO > SrO 
> BaO (lowest). This can be attributed to the fact that 
the polarizing power of the cations increases with the 
decrease of the cation size or ionic radius. The oxygen 
ions bonded to the bigger cations will be highly 
polarized and the replacing of these cations by PbO 
will make the least effect and vice versa. 

When silica is replaced by alumina in these glasses, 
every molecule of alumina will acquire an oxygen from 
the PbO and will form AlO, group, besides the SiO, 
tetrahedra housing the Pb ions in the interstices 
formed by SiO, or AlO, groups. The oxygen ligands 
are not completely screened. Thus, it would be 
expected that the proportion of the Ni ions in the 
octahedral coordination will increase by the 
replacement of silica by alumina. 

When the tetravalent metal oxides, titania or 
zirconia are introduced in these glasses replacing silica, 
the polarizability of the oxygen ligands will decrease. 
Thus, the Ni ion would be expected in the octahedral 
coordination in these glasses. 
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A numerical integration method using Weddle’s rule has been applied to compute the limiting Debye temperature (0) of 47 
cubic and 14 hexagonal crystals from the elastic constant data near OK. These computed values are found to be in very good. ae eed 
agreement with the calorimetric Debye temperatures (5) confirming the equivalence of 0% and 0 in low temperature regions. 
, #2 es re 
lI : : cote spp 
ata 2 Theory -, oye, ak ieee 
_ Because of its simplicity, the Debye model has been For a crystal havi dtnins ee the 
extensively used to explain various thermal, optical — expressi ae ee : ne OS ee sci cell, : sale 
: ae ~vP pression lor Vo can be written as — Ga: ea ears 


and electrical properties of solids. This, in turn, has led : ee ee a FS 
te - 8 : : re h rN4p = ; ity Tes oe Bye Tae J 

© various experimental methods for the estimation of % =7-(9 eae? ak 
the Debye temperature! >. “ Mi aoe ermine 
_ The Debye temperatures derived from single crystal Where h. ky and N4 are the Planck eo 


elast stic constant data near OK (05) acquire special Boltzmann's constant < 
Fe: respectively: p and M denc 
se Sy oes ~{ ; 7) mae s is 7. 
density and the molar weig 


i= [sinoao 


wa J 
UU. gf 


* 


For hexagonal crystals'®, the above expressions can 
be written as 


Ay, =C,,42+0.5(C,, - Cy2)4? + C,,42 

Az, = C45 + 0.5(C,, — C1 2)42 + Cag 22 

Ayy = Oy 342 + Cya(22 + 22) 

Aye =O5(C,, + C,,)2,A, 

Ay, =0.5(C,3 + Cay)Ad- 

Az3 = 0.5(C 13 + Cga)A, A, cachsy 
In the above set of Eqs (4) and (5) C;; (i,j = 1, 2, 3, 4) 

are the elastic constants of the crystals. 


3 Method of Computation 
Because of crystal symmetry, it is sufficient to carry 


a Out the integration in Eq. (2) ina limited portion of the 


- Brillouin zone depending on the structure of the 
crystal. In order that the same computer program can 
be used (with minor modifications) for cubic, 
hexagonal and tetragonal crystals, it is enough to carry 
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out the integration in 1/16 of the Brillouin zone, 
confined in between the directions (100), (1 10) and 
(001). For carrying out the whole computation, a 
computer program in FORTRAN IV was written by 
us. The computer program generates 325 (or 91) 
directions in 1/16 of the BZ and solves the cubic 
equation in 7? derived from Eq. (3) for each direction. 
Then Weddle’s rule is applied for carrying out the 
numerical integration of Eq. (2). The final result is 
printed in terms of the Debye temperature. 


4 Results and Discussion 

Using the computer program, 6 values for 47 cubic 
and 14 hexagonal crystals were computed. The input 
data for the cubic and hexagonal crystals are given in 
Tables | and 2 respectively. In 15 cubic crystals, more 
than one independent sets of data for the elastic 
constants were available. In such cases, computations 
were done for all the available sets of data to see their 
effect on the calculated values of 06. 


Table 1—Input Data for Cubic Crystals 


T.K p,gicm Mg ts (ee Ca 


10'' dyne/em? 


12.296 7.080 3.0920 
11.430 6.1920 3.1620 
11.630 64800 3.0900 
17.400 5600 3.593 
0.259 0.217 0.160 


3.350 1.025 ~—‘1.002 
3.437 1.035 ~—-:0.999 
2737. 80.793 -) 0825+ 
2.725 0.767 
17.620 12.494 
13.150 4.948 


12.107 5.477 
5.660 

16.967 

nou. 

3.788 

8.980 5.025 
8998 $068 


Crystal 1 & K Ps g/cm? M, £ Cia Cis Cus 


10'' dyne/cm? 


39.096 0. 0.341 
74.549 4. 0.540 
58.094 7. 1.350 
58.094 7. ‘1.473 
166.00 3. 0.220 


85.470 0. D288". 4 
165.38 3. 0.474 0.408 
120.92 : 0.676 = 0. 
120.92 4, 0.676 0.493. 
212.37 3.210» 0.360 -~— 0.28 
28.086 16.772 6.498. 
107.87 13. 9.733 
145392 93 pi 
143.32 
22.990 

~ 22.990 

102.90 
58.443 
58.443 
41.982 

0 41.988 

762 149.89 


\ 
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Table 2—Input Data for Hexagonal Crystals 
C13 
10!! dyne/cm? 


0.0 1.8477 18.0244 29.94 2.760 1.100 
0.0 8.8415 224.80 12.923 3.990 4.095 
4.0 8.9332 117.86 31.950 6.610 10.210 
4.0 7.888 314.50 7.400 2.8220 2.034 
4.2 12.773 356.98 19.010 7.450 6.550 
4.2 7.4713 229.64 $392 3.871 4.513 
0.0 1.7752 48.610 6.348 2.594 2.170 
4.0 12.3875 202.14 57.630 18.720 16.730 
4.0 4.5225 95.800 17.610 8.690 6.830 
4.2 11.8170 408.740 4.440 3.760 3.000 
4.2 4.4990 177.82 8.340 2.910 1.900 
4.2 7.2310 130.76 17.909 3.750 5.540 
4.0 6.5280 182.44 15.540 6.720 6.460 

9.9235 349.94 9.2310 3.230 2.800 


Crystal oy p, g/cm? M,g E34 Gea 


Excepting Lu, all the elastic constant data have been taken from 
Simmons and Warg'*. The data for Lu have been taken from Ref. 16. 


Table 3—Calculated Elastic Debye Temperatures (0%) and es 
the Calorimetric Debye Temperatures (05) for Cubic Crystals — 


kK -— “pik Ref. for ‘Crystal, 05, K 
, BPS ina ame Sees 
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Table 4—Calculated Elastic Debye Temperatures (0%) and 
the Calorimetric Debye Temperatures (08) for Hexagonal 


Crystals 
Crystal 0%, K %, K %, Diff. Ref. for 0% 
Be 1463 1440 1.5 20 
Cd 212.5 209 1.6 20 
Co 470.0 345 5.3 20 
Gd 181.3 187 —3.1 41 
Hf 253.1 252 0.43 20 
In 104.8 108 —3.0 20 
Mg 386.1 400 nat 20 
Ru 554.4 600 —8,.2 20 
Ti 425.4 420 1.2 20 
Tl 79.59 78.5 1.3 20 
Y 257.1 280 —89 20 
Zn 327.0 327 0 20 
Zr 260.0 291 1.6 20 
Lu 185.6 196 — 5.6 42 


The results of the computations are presented in 
Tables 3 and 4. In Table 3, the results correspond to 
325-direction integration while in Table 4 the values 
correspond to 9l-direction integration. In these 
Tables, #§ represents our calculated values and 0 the 
literature values. 

On theoretical grounds, it is well known?!’~ '? that 
in the low-temperature limit 05 = 06. It is seen from 
Tables 3 and 4 that in most of the cases 0§ = 06, within 
the limits of experimental uncertainties. The few cases, 
where the differences are seen to be larger than the 
experimental uncertainties, can be explained satisfac- 
torily. For example, let us take the case of Li where the 
difference between 06 and 0 is — 5.7%. It is known that 
Li undergoes a martensitic phase transformation near 
78 K and below this temperature both the bec and hep 
phases exist. The experimental values of 05 correspond 
to a partially transformed phase. Further, we had to 

-_use elastic constant data at 78 K as our input data. 
Thus, we see that a proper comparison between 05 and 
0% is not possible in the case of Li. 

In AgCl, the results of our calculation for 05 for two 
independent sets of data differ by less than 1 %,. But the 
_ difference with the quoted experimental values of 0% is 

13% of 0 which is unusual. This difference may be due 
© an error in the quoted value of 06 given in AIP 
-Handbook?® without any reference to the original 


es comparative large differences for the 

als Co, Ru, Y, Gd, Hf, Mg and Lu point 

ssity of more precise determination of 

cheat as well as the elastic constants at very 
‘atures. On the whole, the agreement 

and Op is bvety 8 ood. 

trova and Fedorov?! have pointed out 
ebye temperature % on the 


C AJ 


ATOMIC NUMBER 


Fig. 1— Variation of 0% with atomic number of elements and 
compounds 


atomic number of the elements in the periodic table. In 
Fig | we have plotted the Debye temperatures of some 
elements and compounds against the atomic numbers. 
It can be seen that there is a regularity in the variation 
of 06 for different elements of a sub-group in the 
periodic table. This dependence of 0% on the atomic 
number may be used for making rough estimates about 
the 0% values of an element or compound. For example. 
it may be assumed that 


Or, < Oe. < Onn < 91K 
where 91 K is the Debye temperature for potassium. 


A comparison of the results of calculation for 05 by 
using the 325-direction integration and 91-direction 
integration showed that the percentage difference 
between the two sets of results was less than the errors 
involved in the determination of the experimental data. 

Our present study leads to the conclusion that the 
numerical method of integration using Weddle’s rule is 
an accurate method for the computation of the Debye 
temperature. Further, it is also concluded that the 
limiting Debye temperatures obtained from elastic 
constant data are equivalent to those dtrived from 
specific heat data. One consequence of the equivalence — 
of 0% and 0% is that the Debye approximation can be 
used for the accurate evaluation of the lattice 

contribution to the specific heat at yey low 
temperatures (7/0) < 0.01). 
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Molecular magnetization, remanent magnetization, Remanence ratio and cation distribution are.computed for the slow- 
cooled and quenched samples of CuFe,O,. A-B interactions are strengthened on quenching the samples. The remanence ratio, 
however, is found to be unaffected by quenching which is explained as due to the impedance to domain-wall motion. All the 
samples exhibit almost the same hysteresis loss. The quenched samples show higher permeability and less coercivity. Quenched 
as well as slow-cooled CuFe,O, samples show tetragonal distortion of the lattice which decreases with the increase of the 


1 Introduction 
CuFe,O, is highly sensitive to heat treatment’ and 
shows semiconductive properties*. Switching and 
memory phenomena are also exhibited by CuFe,O, 
when it is quenched from elevated temperatures’. 
. However, many controversies exist on the behaviour of 
1 CuFe,O, and no conclusive results on any property 
shown by it are available, which renders CuFe,0,4 

interesting for further studies. 
The cation distribution in a ferrite is usually found 
from various studies, e.g. X-ray diffraction‘, 
Mossbauer effect®, and magnetization. However, 
studies on the sensitivity of cation distribution and 
_ magnetization to heat treatment have given rise to 
‘ > some controversies and inconclusive results. These 
difficulties can be overcome from studies on 
magnetization of samples taken through various bose 
cS I atme nts. 

= inthis paper, we report our work on magnetization 
on slow-cooled and quenched samples of 


ing 3 temperature on tetragonality, which is 


i 


connected with Cu?* Jahn-Teller ions on B- 


ne CuFe,O, was prepared using AR 
. and Fe,03 by the standard ceramic 
pletion of the solid-state reaction was 
ne help. of X-ray diffraction studies. 


: 950°C for 4 hr and cooled i in 


Q quenching oft the remaining 


er and 2 mm thickness were > 


' temperature of quenching. This is explained by cation transfer. 


High-field loop tracer HS 869, supplied by 
Electronics Corporation of India, was used for the 
measurement of magnetic saturation (M,), remanent 
magnetization (M,) and remanence ratio (M,/M,). 
Measurements were carried out on CRO screen which 
was properly illuminated and calibrated. M, at various 
elevated temperatures was computed using Hs 868 A 
hysteresis and susceptibility apparatus. The samples 
were prepared in the toroid form for this measurement. 
Curie temperatures for these samples were determined 
by the technique described elsewhere’®. 


3 Results and Discussion 

Fig. 1 shows the variation of molecular 
magnetization (ng) with temperature of quenching for 
CuFe,O, samples. It is observed that ng is maximum 
for the samples quenched from 800°C and minimum 


@ Slow cooled sample 


Molecular magnetization (ng), Bohr magneton 


500 600.700 800 
a Temperature of quenching, °C 


Fig. 1—Variation of ng with temperature of quenching 
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for those quenched from 600°C. The point marked by 
solid dot shows value of ng for the slow-cooled sample. 

The values of ng, M,, M,/M, and lattice parameter 
ratio (c/a), energy (E) and hysteresis loss at room 
temperature are given in Table 1; the values of M,/M, 
have been calculated from studies on hysteresis and 
those of c/a from X-ray diffraction studies. The values 
of cation distribution, Curie temperature, g-factor at 
OK are given in Table 2. Values of ng at 0K have been 
calculated from M, versus T plot (Fig. 2) by the 
extrapolation method. 

The value of ng for slow cooled sample is in close 
agreement with the earlier reported value of ng for 
CuFe,O, (Ref. 7). The minor deviation may be 
attributed to the structural sensitivity on account of 
factors like porosity, defects, conditions and 
atmosphere of firing and the attendant changes in the 
cation distribution. 

Values of g-factor are calculated from M, versus T 
plot®. It can be-concluded from these values of g that 
the possibility of Cu* being present in Cu-ferrite is 
almost negligible and so the cation distribution given is 
justifiable. It is apparent from Table 2 that as the 
temperature of quenching increases, there occurs more 
migration of Cu** ions from B-site to A-site. As a 


Table 1—Values of ng, M,, M,/M,, c/a, E and Hysteresis 
Loss at Room Temperature for Quenched and Slow Cooled 


CuFe,0, 
Thermal Np M, M,/M, c/a iE Hyste- 
history (Bohr (emu/g) eV __ resis loss 
magneton) (Arbi- 
trary 
units) 
Quenched 1.5 19 0.54 1.037 0.044 215 
from 800°C 
owe 1.40 18 0.55 1.039 0.052 221 
trom 700° 
Songrad 120 16 = 057 1.042 0068 94 “200-100 0 100 200 300 400 
=~ cooled ‘1.10 14 0.56 = 1.052 0.094 9 Temperature.°C 
9 c . . . 
rom 950°C Fig. 2—Variation of M, with temperature of the sample 
Table 2—V istributi 
alues of ng (at OK), Cation Distribution, g-factor and Curie Temperature for Quenched and Slow Cooled 
CuFe,0, 
Thermal Np at OK , ee 
history (Bohr Cation distribution g-Factor Curie 
magneton) — 
Slow + ' CO) 
pee 1.20 (Cus .o25 Feo 575)4 (Cua 575 Fe? $,5)"O, 2.13 440 
Quenched 1.50 2+ - 
ton crc (CUd 0615 Fea.5395)4(Cud bags Fe? b61s)°O, 2.13 415 
Quenched : + - 
irom 700°C 1.80 (Cus 4 Fe} 3)" (Cug.5 Fe} +)" O, 213 390 
Quenched 2.00 2+ : 
from 800°C (Cus 25 Fes 875)" (Cua },, Fe} 425)"O, 2.13 360 
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result of this the Curie temperatures are lowered and 
A-B interactions are strengthened. Therefore ng shows 
an increase. 

It is seen from Table | that M, also increases with the 
temperature of quenching. Locking of domain-wall 
motion due to pores and other hindrances may be 
responsible for this behaviour of M,. The change in M, 
also indicates that the temperature of quenching has a 
bearing upon the impediments to the domain-wall 
motion, although this change in M, is very small. 

M,/M, has not been found to vary with quenching 
temperature within experimental errors. 

For CuFe,0O,4, higher the quenching temperature 
higher are the values of M, and M,, yet M,/M, depends 
on factors like impurities, defects, polarizable 
constituents and other hindrances to domain-wall 
motion which determine the retentivity. The 
hindrances to the domain-wall motion are annealed 
out at 800°C to their equilibrium concentration and so 


quenched fro 


a 


a 8 
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| 


ng, Bohr magneton 


+0 


4035 1037 1039 1042. 
s Axial ratio c/a 
| Fig. 3—ng versus c/a 


+045 


it in less degree [as is evidenced by the 
of the coercive field (H,)] in the samples 
1 from 800°C because the defect clusters are 
avoured at higher temperatures; annealing out of 

sf ences to the domain-wall motion at 600°C is 
of defect cluster formation at lower 
| — so the sample quenched from 600°C 
NV ss tnt a hence M,/M, does not show 
ice on the ching temperature. 
g ne piece ile ions from B- 


sn eee iw am. | 
= a ae a arr ; 


permeability of the samples increases while the 
coercivity decreases. The saturation field for all these 
samples is found to be less than 20 Oe. 

Fig. 3 shows the variation of ng with axial ratio (c/a). 
It is seen that as the tetragonal distortion decreases, ng 
increases almost in linear proportion. This is because of 
the change in cation distribution on quenching, as 
discussed already. 

The samples of CuFe,Ox4, slow cooled from 950°C 
and air quenched from 800, 700 and 600°C, showed a 
tetragonal distortion of the lattice. It is seen that, as the 
temperature of quenching changes, the tetragonality 
also changes its degree. The axial ratio c/a is maximum 
for the sample quenched from 600°C while it is 
minimum for the sample quenched from 800°C. In 
literature, c/a for CuFe,O, is reported to have values 
from 1.03 to 1.06 (Refs 9-11). Our values of c/a show 
close agreement with those reported in the literature. 
The tetragonal distortion of the spinel structure is the 
result of sufficient concentration of distorting Cu?* 
Jahn-Teller ions present on B-site, as evidenced from 
Table 2. On quenching, the cation distribution is 
altered and frozen in, showing variation in axial ratio. 
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Sharan and Prakash [Indian J Pure & Appl Phys, 5 (1967) 442] from their study pea halides — gee ee 
aration at the point of inflection to that at absolute z 
Mayer’s potentials concluded that the ratio of the interatomic separa ; 
aN constant. The same has been used for getting the ratio Up/kT,, (where Uo is the interaction energy per particle at 
absolute zero of temperature and T,, the melting temperature) and also the Sutherland-Lindemann empirical melting point 
equations from Mie’s potential. It has been concluded that the Sutherland-Lindemann constant is structure dependent. The 
agreement between the theory and experiment for Uo/kT,, is good for alkali halides. The theory has been extended to include 


order-disorder effect in binary systems. 


1 Introduction 

The equation for the average potential energy of a 
particle with its equipartitional thermal energy near 
the melting point gives the famous Sutherland- 
Lindemann equation. Mathematically it can be 
expressed in the form 


Vo=LV 213 7 1/2 y-1/2 ET) 


where Vo is the frequency of vibration of the particle 
and is often taken to be restrahlen frequency; L, V,,, Tn 
and M are respectively the Lindemann constant, molar 
volume at the melting temperature, melting 
temperature and the molar weight of the particle. By 
substituting experimental values, L is found to have an 
approximate value of 2.8 x 10'? cmg!!? sec”! deg~ 1/2 
The value of the constant has been fixed completely in 
an empirical manner. A critical analysis would show 
that strictly speaking L should be structure dependent 
and differ from substance to substance and, therefore, a 
theoretical derivation which would include the above 
variations is called for. Dheer and Sharan’ from their 
analysis on the isotherms of caesium bromide pointed 
out that the moment the interparticle separation ina 
binary system exceeds the separation corresponding to 
the point of inflection of the potential energy versus 
interparticle distance curve, the crystal structure would 
collapse. The criterion has been used for alkali halides 
and metals by Sharan and Prakash? for the 
determination of the ratio of the interparticle distances 
at the melting point to that at absolute zero. The 
dependence of the Product Po (P and » are pressure 
and volume respectively of the substance) on 
temperature, internal energy, which in turn can be 
expressed as a function of frequency of vibration 
coupled with the melting criterion offers a possibility of 
theoretically deducing the Sutherland-Lindemann 
equation. The derivation of this equation for a 
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monatomic solid obeying Mie’s potential and its 
application to alkali halides forms the subject matter of 
this paper. 


2 Derivation of Sutherland-Lindemann Equation 
According to Mie®, using the concept of second 
virial, the equation of state can be written in the form 


1} dU 
ie aicee see 


where U is the interaction energy per particle. The 
average value of the third term in Eq. (2.1) for Mie’s 
potential energy 


ee 


5 ee 
stall 
coupled with the equipartition law and the equilibrium 
condition r= ro, reduces Eq. (2.1) to the form 


Py= Fn +m+1]}+ =a — {@y — eae 
3 |n—m r r 
aR fee 
where the symbols have their usual meanings and the 
subscript zero refers to the value of corresponding 
quantity at the absolute zero of temperature. The ratio 


of r,, and ro, the interparticle separations at the melting 
temperature and at absolute zero respectively, can be 


expressed as 
n+1 }4. 
. Es i ...(24) 


where r; is the interparticle separation at the point of 
inflection in the U versus r Curve and is given by 
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Sr -,, 70: The constant K can be determined 


according to the criterion develo Sharan 
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Fig. 1. Writing Eq. (2.3) for melting temperature T,, PF. 
and v,, and substituting r,,/ro from Eq. (2.4), we get 
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The constant L depends upon the nature of the 
substance and its structure through o, mand nand can 


be theoretically evaluated by the formula 
1 


2 


kKN*3[n+m+1]{n-—m] 
| — — 4:3 


m+2 n+2 
an'| l feet l ia 
K"™*?*I n4+1 K**41 434 


3 Results and Discussion 

The value of L for potassium chloride using K 
= 0.8625, n= 9.0, m= 1.0 and o =2 comes out to be 
2.641 x 10'* cm g!’? sec” 'deg™ 2. 

The value of L Pest, in this manner is of the 
right order of magnitude and is about 5 °% lower than 
the accepted value of 2.8 x 10'? cm g!/?sec™ !deg~ 1/2. 
Though ZL is structure- and substance-dependent, the 
value is not significantly different from the accepted 
one. As a further check, Uo/kT,, has been evaluated 


...(2.10) 


Table 1—Theoretical and Experimental Values of — Up/kT, 


Alkali halide — U,/kT,, 

Theor. value Exptl 

by the present 

author 

LiF . S230 — $5.34 aes 
LiCl 39.75 57.56 
LiBr > 40.35 58.76 
Lil 40.75 62.69 
NaF 40.47 opi, 
NaCl 250 Sea 
NaBr A 42.81 ; a 3:70- 
KF 42451 1) aia ie 
KCl 457) Saba 


* gt 
% tS “ne 
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from Eq. (2.7) by assuming I-n potential for the alkali 
halides. The agreement with the experimental values iS 
good except for lithium halides (Table 1). 

The dependence of L on the power n for different 
structures, viz. fec, bec, NaCl and CsCl, is shown in 
Table 2. The consideration of disorder effect in the 


lattice changes L slightly. 


4 Contribution to L Due to Disorder Effect at 7, 
It has been derived theoretically as follows. 
According to the theory of order-disorder, the free 
energy component As — As=;, for an equimolar binary 
system, which we take to be formed of atoms and holes, 
is given by 
Ag — As=; = A= NkT[(1 + S)In (1 + S) 
+ (1 —S)In(1 — S) —2In 2] 
ee +4No(1 —S?) .. (4.1) 
where S is the degree of order (S = 1 corresponds to the 


completely ordered state) and ¢ the interaction energy 


of disorder. For an isothermal system containing a 


a 


constant amount of pure substance, T and N are fixed 


Fa 


— nd Ais a function of S and ¢ only. Then 


E f eel  4¢\- 
pee -3). 


ting this in Eq. (4.2), 
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- where Qo and & are constants. 
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A thermodynamic perturbation method based on Gibbs-Bogoliubov inequality has been applied to a non-simple liquid 
metal, viz. Bi through Heine-Abarenkov-Animalu pseudopotential and this procedure yields the possibility of two packing 
densities. Static structure-factor calculation with two packings yields the shoulder characteristics quite reasonably. 
Calculations at the close packing of the liquid yield satisfactory results. 


‘1 Introduction 

Successful calculations’ ~ ° of the thermodynamics of 
simple liquid metals and their alloys have been possible 
by the combination of the pseudopotential theory and 
thermodynamic perturbation theories through a 
variational method based on the Gibbs-Bogoliubov 
(G-B) inequality®. Applications of these methods to 
non-simple liquid metals like Bi have been considered 
to be not possible’. In this paper, we have however, 
shown that the pair potential description is not 
unsuitable to liquid Bi and this procedure can account 
for some peculiar properties of this liquid. The Heine- 
Abarenkov-Animalu (HAA) model pseudopotential® 
together with Vasistha and Singwi (V-S)’ exchange and 
correlation function provide necessary energy 
wavenumber characteristic and the corresponding 
effective pair potential. A reference system of hard 
spheres in the Percus-Yevick (P-Y) approximation’? is 
used, the hard core diameter being chosen through 
usual procedure. This ab initio calculation yields 
several interesting results which are discussed in Sec. 3. 


2 Theory and Calculation 
In the pseudopotential approach, the total energy of 
a metal may be expressed as the sum of a volume- 
2 dependent term independent of ion positions and other 
terms which vary with structure. In terms of 
1ormalized energy wavenumber characteristic Fy(q), 


t e effective pair potential energy g(r) and pair 


ae | . Fu(q) nt dg .() 


(1 4\ bw? 
a) i = f(a) 
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where Z* is the effective valency, Q the atomic volume 
and f(q) accounts for exchange and correlation effects 
of conduction electrons’. The term €(q) is the dielectric 
function of interacting electron gas'’. The term w°(q) is 
the HAA potential in g-space®. 

The free-energy may be expressed in terms of the 
packing density 7 of hard spheres as detailed in the 
following. 

F,,=3keT —TS,, ee (2 
The first term being internal energy of an ideal gas and 
in the second term the entropy S,, is given by, 

Shs = Sgas +S, Yet 
Sgas being that of the gas and S, is the packing density- 
dependent part and is obtainable from the reference 
system. According to G-B principle, the free-energy for 
the system, 

Pa PF, +<M%,) ive (4) 
<V,;> being the configuration average of the total 
interaction potential performed over the reference 
system. 

Since the inequality (4) constitutes the upper bound, 
it follows then by treating 4 as the variational 
parameter that the best hardcore diameter can be 
chosen by putting 


OF 
(Fan ® ..(5) 


Another important consequence of Eq. (5) is that 
S= Shs ss Sec 


2 
_ where Sace = 7 N(Es)k5T 


and 
Su= Sys + Nia tn —n)+ 34! is eal 
Also, 
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3 Nkp X 2n(2— + 7/2) ( Olnn 
Ca=5Nka—-—— "j= din T Je 
+ Coelec 


Ca is the specific heat at constant volume and other 
terms have their usual meaning. C,, the specific heat at 
constant pressure can be obtained by numerical 
differentiation of the tota! entropy. 


3 Results and Discussion 

The well known variational procedure yields two 
possible hardcore diameters satisfying Eq. (5). At 
higher temperatures consistent with the usual 
expectation the two diameters gradually tend to 
approach one another (Fig. 1). This is really interesting 
and is not known for simple liquid metals. At a given 
temperature and volume, two configurations are thus 
possible. The two collision diameters are not unlikely 
for these types of liquid metals'? which show shoulder 
on the right hand side of the structure factor peak. In 


awe> For.é"s: 


T,K 


Fig. 1—Variation of n and o with temperature [----o’s; 


n’s] 


Wir) 210% au. 


P (r) x10 a.u. 
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Fig. 2—Pair potential of liquid Bi at m.p. (Inset: Pair interaction of 
liquid Bi at m.p.] 
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Fig. 3—Structure factor of liquid Bi at 271.3°C 


Fig. 2 the effective interionic potential deduced from 
Eq. (1) has been shown. In Fig. 2 [inset] we have also 
shown the corresponding force curve. The force curve 
evidently shows two possible collision diameters. It is 
expected that at low temperature, the core diameter is 
indicated by the point A while at higher temperature 
indicated by the point B. Average thermal energy at 
melting temperature corresponds to the point C. Heine 
and Weaire’* have shown that the Coulomb forces 
favour close packing and that other cases may be 
explained by special circumstances. Thus if the 
nearest-neighbour distance for close packing occurs 
near the energetically unfavourable point B, it helps 
some such neighbours to move inwards while others to 
move out so as to conform to energy requirements. 
Heine and Weaire’* suggest that similar con- 
siderations have an important role to play in liquid- 
State structure theory as well, specially in the case of 
those metals whose volume diminishes on melting. 

A calculation by Silbert and Young!‘ is also relevant 
in this connection. They performed a model calculation 
for hard spheres with repulsive tail (like two hardcore 
diameters) and could reproduce very nicely the 
Shoulder of the structure factor in liquid Bi. 

Orton'® has suggested a model for such liquid 
metals based on two near neighbour distances. Thus Bi 
molecules in the liquid may come together to give _ 
either long or short neighbour distances. The majority — 
of the molecules have the longer near neighbour 


distances, while a minority of the lecules. 
short distances. On the henia otis model ae sould - 
calculate the structure factor and could produce the 
shoulders in this type of liquid metals. We have 
calculated the structure factor on thé basis of Orton’s 
model with the two packing densities obtained from 
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Table 1—Variation of Sound Velocity with Temperature 


— Packing density Sound velocity (c), m/sec 
"a "es Calc. Obs. 
544.3 * 0.483 0.429 1577.8 1647.0 
573 0.480 0.427 1566.7 1618.8 
773 0.460 0.423 1574.0 1649.1 
823 0.450 0.420 1566.9 1639.7 


in the long wave limit has been obtained on Orton’s 
model in the hard sphere approximation and this hard- 
sphere value has been used in the calculation of sound 
velocity of liquid Bi based on Ascarelli’s method'®. The 
constant in the Ascarelli’s expression has been 
obtained through compressibility sum rule. At higher 
temperatures, we have the required packing densities 
through the variational method. The atomic 
concentration was assumed to be proportional to 
packing density (core diameter was assumed to be 
constant in the small range considered). The results 

(Table 1) show that the sound velocity in Bi remains 

almost constant with temperature. However, they 

indicate a tendency to decrease with increase of 
temperature which is in conformity with experiment. 

Experimentally'’, temperature coefficient of sound 

velocity is negative and almost zero for Bi near the 

melting temperature. In view of the limitations of the 
method one cannot expect better agreement. 

We have calculated the excess entropy for liquid Bi 
for both types of packing (— 3.546 for n = 0.429 and 
—4.656 for n =0.483). The experimental value lies 
midway between the two calculated values. The 
= s specific heats at constant volume and at constant 
ao: _ pressure for close packing of this liquid are also 
= Bepemable (Table 2) and their ratio (y) (Table 2) is in 
| ‘good agreement with the experiment. We have 
- rs ‘ > obtained the isothermal compressibility (Table 2) 
< for larger packing from the free energy. 

: reewine Young!®, we next calculated the softness 

_ of the core corresponding to close packing of the liquid. 
i AlnQ 

ainT ), and WaT hse are 0.041 and 

‘res spectively. e = difference between these values 

‘s hows that packing increases considerably on 

7 cted for these types of liquids. 

9 as estion is why there be two such packings 
3i7 The reason may be that the Cauchy 
rated by the electron gas is sufficiently 

ions close to each other in liquid state 

1ts of the low pressure solid phase- 

! in some parts of the liquid. 

there i is a partial retention in 

valent or homopolar 
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Table 2—-Thermodynamics of Liquid Bi at Melting Point at 
the Close Packing (yn = 0.483) 
Parameter Obs. Calc. 
Excess entropy — 3.90 — 4.656 
(Nkp) * 
Isothermal 4.14 2.424* 
compressibility 
(B,) x 107 '*cm?/dyne 
C,, cal/g 0.0316 0.0209 
C,, cal/g 0.0363 0.0241 
C,/C, =y 1.15 1.15 
Vp 2.00 2.67 


*From free-energy 


in a metallically bonded matrix. Further work on this 
point is needed before one can conclude anything. But 
what is important is that two possible packings are 
likely in a liquid like Bi where a decrease in volume 
occurs On melting and consequent increase in 
electronic pressure. As a further corroboration of our 
results, we calculated the pseudo Griineisen constant 


- z. oa using the experimental thermal expansion7°® 
IT 

and other calculated quantities. The value so obtained 

is 2.67 while the experimental 7, comes out to be 2.00. 

Using the theoretical Griineisen constant of 2.67 and 


the Debye equation, 


Feta 5s 
c= | Gp? | = 21,45 


we get C, = 7.0,a reasonable value for metallic liquids. 
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The velocity of ultrasonic waves (2 MHz) for different compositions of mixtures of CCl, + o-toluidine and CCl, + p- 
toluidine has been measured in the temperature range 10-50°C and 30-55°C respectively. The velocity (c) versus composition 
(C,,) plots exhibit peak at ¥,: ¥, =8:2in CCl, + p-toluidine. The excess compressibility BF, excess van der Waals’ parameter 
b¥ and(dc/dT)* have been diseussed in the light of intermolecular AB interaction and resulting disorder in these mixtures. It is 
observed that the strength of AB interaction depends not only upon the nature of the groups but also on their relative positions 


in the ring. 


1 Introduction 
Ultrasonic parameters are being extensively used to 
study molecular interactions in pure liquids’, binary 
and ternary liquid mixtures” ~°, and ionic interactions 
in single and mixed salt solutions’’*. A departure from 
linearity in the velocity versus composition behaviour 
in liquid mixtures is taken as an indication of the 
existence of interaction between the different 
species> ~*. However, it is seen that a representation in 
terms of the observed parameters, such as velocity of 
the ultrasonic waves, has a limited utility. Such a 
representation does not provide any information 
about the nature and the relative strengths of the 
various intermolecular interactions. On the other 
hand, a number of theoretical’ and experimental'° ~'? 
investigations have shown that a representation in 
terms of the derived parameters, such as adiabatic 
compressibility (8), van der Waals’ parameter (b) and 
their deviations, such as B®, b® (excess compressibility, 
excess b, etc.) from those given by additive rule, 
provide a better insight into the intermolecular 
processes. Interactions in binary liquid mixtures have 
been studied by various authors'*~'’ from the 
knowledge of excess thermodynamic functions and 
acoustical parameters. It appears that >a many 
attempts*''® have been made to use ultrasonic 
: dae to test the liquid state theories’? pin 
very few investigations have been made!?:!2-22:29 to 
determine relative strengths of homo- and heteromole- 
r interactions in binary liquid mixtures. This paper 
ents the results of an ultrasonic study of mixtures 
7 tetrachloride with ortho- and para- 
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2 Experimental Details 

The liquids used were BDH Analar grade and 
distilled thrice before use, the middle fraction being 
taken in each case. The mixtures were prepared 
immediately before use, by mixing appropriate 
volumes of the constituent liquids taken correct to 0.1 
ml, the total volume being about 100 ml. The 
ultrasonic velocity was measured at 2 MHz by the 
interferometer method using Mittal’s M’77 instrument 
in the temperature range 10-50°C and 30-55°C for 
CCl, + o-toluidine and CCl, + _ p-toluidine 
respectively. The details of experimental technique 
were the same as reported earlier’’. 

The parameters B, and b were obtained using 
standard relations?*, and (dc/dT) was obtained from 
velocity versus temperature plots. The excess 
parameters BE, b= and (dc/dT)* are given by the 
difference between the observed and theoretical 
(simple additive rule) values of the respective 
parameters. The error for each of the excess parameter 
was estimated and the deviations from the additive law 
were found to be outside the limits of error (* 3%). 


3 Results and Discussion 

The velocity and compressibility versus composition 
plots are shown in Fig. 1. The f¢, b*; V* and (de/dT)* 
variations are presented in Figs 2-5 respectively. The 
variation of density versus composition in both the 
systems is nonlinear. 

The BE variation (Fig. 2), in CCl, + o-toluidine, and 
CCl, + p-toluidine mixtures is negative over the whole 
concentration range and indicates the presence of an 
AB interaction. A large dip in Bf at ratio 8:2 is 
observed in CCl, + p-toluidine, which reflects strong 
AB interaction leading to complexation. The peak in 
ultrasonic velocity and corresponding dip in density 
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Fig. 1—Velocity (c) and compressibility (B,) versus mole fraction 
(Cm) of o-and Pt at 40°C [Curve: 1,3, o-toluidine; 2,4, p- 
toluidine] 
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Fig. 6—Density (d) versus mole fraction (C,,) of p-toluidine.at 40°C 
the ring. This increased dipolar character of the isomer 


interaction among the meta- and para-isomer. The 


tendency of the para-isomer to form complex is due to 


the = cape between the oe in pres ng 


is responsible for the observed larger strength of AB 
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toluene. It has been found? that NH, group 
interaction with CCl, decreases the effective molecular 
size, while CH, group interaction with CCl, increases 
the molecular size. The large peak in b* for CCl, + p- 
toluidine at the ratio 8:2, and b* variation for CCl, + 
o-toluidine beyond the ratio 8: 2, reflect the major role 
played by CH, group in AB interaction. This 
interaction of CH, group seems to be maximum at the 
ratio 4: 6, evident from the large peak in b* observed in 
CCl, + toluene, CCl, + m-toluidine and CCl, + 0- 
toluidine mixtures. Such a peak is not observed in CCl, 
+ p-toluidine mixture. This may be due the balance of 
CH, and NH, group interaction. The overall negative 
b® variation in CCl, + toluidines also shows the 
important role played by NH, group in AB 
interaction. Thus both CH; and NH, groups take part 
in AB interaction in this system. 

The prominent feature observed in all the systems is 
the close resemblance of b* and V* parameters (Figs 3 
and 4). It is worthwhile to mention here that V™ in 
systems CCl, + aniline and CCl, + o0-toluidine is 
negative over the whole range of concentration. 
Adgaonkar et a/.*® have reported positive V® in the 
same system. The present negative b* variation over 
the whole concentration range, however, corroborates 
the V* systematics in our study. Thus it is apparent 
that the evaluation of b lends support to the 
conclusions drawn from molar volume (V) 
systematics. 


4 Conclusions 

The present study on CCl, + o-toluidine and CCl, 
+ p-toluidine shows that (a) the AB and BB 
interactions are of nearly equal strength, while AA 
interaction is very weak. (b) The AB interaction, in the 
nature of hydrogen bonding, leads to the formation of 
complex molecule in CCl, + p-toluidine of the type 
4CCl,-C,HoN. (c) The strength of AB interaction 
depends not only on the presence of groups but also 


upon their relative orientation in the ring, decreasing 
in the order 0, m, p. 
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Application of HHOB Approximation to e~-Li Elastic Scattering 
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High energy higher-order Born (HHOB) approximation as proposed by AC Yates [Phys Rev A (USA), 19 (1979) 1550] is 
applied to calculate the elastic differential cross-sections and total collisional cross-sections for e~-Li elastic scattering at 
intermediate energies. Results of the calculations are compared with recent theoretical results of F W Byron and C J Joachain 
[Phys Rev A (USA), 4 (1981) 1817] and the experimental data of Williams et al. [J Phys B (GB), 9 (1976) 1529, 1576]. 


1 Introduction 

The HHOB approximation proposed by Yates! is 
one of the successfully applied high energy, small angle 
approximation. This HHOB approximation yields 
reliable results for elastic scattering processes? ~*. In 
the case of e -He atom scattering, results? obtained 
using this approximation are very encouraging. 

In the present paper, we have calculated the elastic 
differential scattering cross-section (DCS) and total 
collisional cross-section (TCS) for e -Li elastic 
scattering. The first and second terms in the Born 
approximation are calculated within the framework of 
HHOB approximation’, the third term is calculated 
using the Glauber eikonal series (GES) of Yates*. We 


as suggested by Walters®. It was investigated by 

Mathur et al.’ that at both intermediate and high 

energies, the differences between the scattering 

parameters obtained by using one- or 'three-electron 

wavefunctions of Li atom were not very appreciable. 

Keeping this in mind, we have used the one-electron 
_ wavefunction for the Li atom. 


2 Theory 
Atomic units are used in this study. K;, K, and q = K; 
_—K,, represent the incident, final momenta of 
scattered electron, and the momentum transfered to 
the target atom, respectively, during the collision 
process. The DCS for fixed q through O(1/k?) can be 


boi [FM +|Im F)|2 + |Re, F2)/2 
 +2F (Re, F? + Re, F? + Fes] seh) 
here F® is the first Born approximation, Im F, 
_ F?) and Re, F are the imaginary and real parts 
er (1/k;) and order (1/k?) amplitudes in the second 
sroximation respectively, and F};, the third 


sed 2) 


have used the one-electron wavefunction for Li atom > 


is derived using Yates°. The matrix element — 


where V is the interaction between the incident electron 
and the Li atom and is given as: 


1 1 
Vo(ro, r)= —— + 


—— + V, Bee (| 
ro WWo-F,| (3) 


where 


V.= | 2 + 27 lexp(- 5.419) 
0 


V. is the core potential and ro, r,, are the position 
vectors of the incident and target electrons, with 
respect to the target nuclei. The wavefunction used for 
the ground state of lithium atom is 

Y,,(r) = Arexp(— y,r) + Bexp(— y3r) ...(4) 
where A = 0.11252; B= — 0.42204; y, = 0.65; y, =2.7. 
In the calculation of higher order terms, we have 
neglected the core potential V, contribution. It was 
shown by Guha and Ghosh* that V, will not affect the 
cross-section appreciably. The closed form of the 
amplitude factors in the HHOB approximation are 
given as: 


2 = 1 
(1) — A =o 
," p> ?» 152 @ +92)” SG? +9) 
nee ee 1 | 
ya(q? + ya)? Yalq? + ya)? (5) 
3 B 1 2 
a. ae Pn — 1° 
7 p> E Pat (q +y2)! y] ..(6) 
3B 1 q? . 
7 Bat | -+—3 18 
Re; r 2 E Dava {21 (q + Yn) s}] se .(7) 


R Pa=> ee ee +51} 
van Lake “OB +2) fd) 


Fls= — 31 arent ee( Sth 
+5 -24@ 99} | 0) 
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where : : : 
q T = q 2 
2)—=2(log—} +— + —-]}]/m 
A(q, Yn) 2( og 4) 6 py ( «) 


Yn 
a) 2\m B 
=-y (23) |e when — > 1 
m=1\4 Ls 

where me 

ep FE _ OO 5, 
2 ee eae? eg 
: fe 10 

where AE is average excitation energy obtained from : 

=) 


Vanderpoorten”. 

A's, Bis, Cis and y,s are constants with following 
values : 

Ap = 8.95318 A, =4.7740 A,=0.63640 _ 0 
B, =0.71247 B,=0.37990 B, =0.05064 

C, = 29.73803° C,=2.49206 C,;=0.00710 

we oe py, = 34 74 =13 

Wi=54 yp =34 y3=13 yy=54 es 
The typical integrals I,, I? and 1,, 19, 1;, I$ are 


analogous to Yates‘. The total cross-section can be etre “Scar rennet] 
calculated using the optical theorem: | Fig. |—DCS for e “Li elastic fer 20 


ee (q=0) se (10) 


sal theorem, we have ‘ie Gatetinead the — 
cross-sections for elastic scattering of = "E pe 


7 
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calculations in HHOB approximation with and 
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i seahy ion nand eikonal Born series results taken from Ref, 
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Fig. 4—DCS for e -Li elastic scattering at 200 eV [The —— is 
same as for Fig. 3.] 
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The far infrared (20-500 cm~'), infrared (200-4000 cm~') and laser Raman (100-3500 cm~!) spectra of 2,5- and 2,6- 
dimethylanilines have been recorded using Fourier far IR spectrometer, Perkin Elmer grating spectrophotometer and He-Ne 
laser Ramanor spectrometer respectively. The vibrational spectra have been analyzed assuming C, and C,, point groups 
respectively for these two compounds. The assignments for the fundamental vibrations, combination and overtone frequencits 
and the internal modes of vibrations of methyl group as well as amino group have been proposed. The vibration number in 
Wilson's notations for each mode of fundamental vibration has been indicated. 


1 Introduction 
The infrared spectra of 2,3- 2,4- and 3,4- 
| dimethylanilines have been studied by Prasad! in the 
region 400-4000 cm '. The infrared spectra of 2,5- 
) dimethylaniline have been reported by Varsanyi? in 
the region 550-4000 cm ' which covers only 70%, of 
the expected normal modes. However, the Raman 
spectra of the same compound have not been 
investigated so far. The vibrational study of 2,6- 
dimethylaniline has not been reported so far. The aim 
of this work is to obtain all the vibrational frequencies 
E in the region 20-4000 cm‘ and to propose complete 
___ assignments for these two molecules. Hence the FIR, 
__IRand laser Raman spectra are analyzed for these two 
molecules and the results are reported in this paper. 


2 «Experir fiend ure 
oy Saige 2,5- and 2,6-dimethylanilines (in 

i € at room temperature) were obtained from 
Laboratories Ltd, England. These 
ed pure samples were used as such without 
r purification. The FIR ep were 
orc led on Fourier far IR spectrometer (model 
"AR 30): [Geseatetorante . The spectra 


idk ied on Perkin E int som 
ere eae otometer in the 


are Teale 


using He-Ne laser beam. The 6328 A line was used to 
excite the spectra. The accuracy of measurements was 
estimated to be within +5 cm‘. The samples were 
held in a glass capillary tube and an interference filter 
was used. The spectra were recorded with incident 
beam parallel and perpendicular to the plane of 
polarization separately and the depolarization ratios 
were calculated for the eianbeecey lying in the region 
100-1800 cm '. 


3 Results and Discussion 
According to Randle and Whiffen*, a group of atoms — 

acting as a substituent (CH; and NH, in the present 

case) may be treated as single mass points and the total 

number of fundamental vibrations will remain 30, 

besides the internal modes of vibrations of the 
substituents and combination and overtone frequen- = 
cies. Assuming the molecule 2,5-dimethylaniline as 
planar, it would belong to C, point group. Out of the 30 
normal modes of vibrations, 21 will be planar (a’)and9 
will be nonplanar (a’). The frequencies of the observed 
fundamental bands, aera ome assil 
along with the vibration number in V Jilson sn ot 
are given in Table 1. The ara 
bands for the same 
along with their nal 
96 CH and NH groups ar ed in a 


if " 
, f 
Sa . 7 
) Tl 


| pecies a4; a , by and b aber 
sae ye 8 


paphaniar. ers earth 4 


- aa ° wae eee an 


IR 


a’ 
vj 3080(4.0) 
v, 3048(5.0) 
v3 3005(7.5) 


vs 1578(8.2) 
vg 1510(8.2) 
vz 1460(8.2) 
vg 1310(7.7) 
Vg 1300(7.5) 
Vio 1220(6.0) 


¥,,; 1168(5.5) 
viz 1150(6.0) 
v3 1085(4.5) 
V4 240(5.0) 
Vis 762(6.0) 
Yi6 565(3.8) 


Ver 490(3.6) 
Vig 448(4.8) 
Yio 290(3.8) 
Vyo 215(3.1) 
Vz, 185(3.4) 


a’ 
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Table 1—Fundamental Frequencies (in cm 


Raman Vibration 
number 
2,5-DMA 
3082(w) 2 
3040(w) 20, 
3010(w) 20, 
1638(s) 5 
1588(m) 8, 
1500(w) 19, 
1460(m) 19, 
= T 
1300(s) 3 
1224(s) 14 
1170(s)p 13 
1150(w) 15 
1092(m) —:18, 
943(s)p 4 
763(s)p 12 
575(s)p 1 
493(s)p 6, 
450(s)p 6, 
290(s) 9, 
215(m) 18, 
9, 
805(s\dp I 
733(s)\dp 4 
550(s)d 16, 
528(s\dp 16, 
-375(s)dp ‘10, 
230m) 5 
ines 17, 
is 


Assignment 


v(C—H) 
v(C—H) 
v(C—H) 
v(C—C) 
v(C—C) 
v(C—C) 
v(C—C) 
v(C—CH;) 
6(C—H) 
v(C—C) 
Kekule vibration 
v(C—NH,) 
6(C—H) 
6(C—H) 
v(C—CH;) 
6(C—C—C) 
v(C—C) 
Ring breathing 

6(C—C—C) 
6(C—C—C) 
6(C—CH;) 
6(C—CH;) 
6(C—NH,) 


»(C—H) 
ECE) 
WE—C=0) 
"WC=C—CF 
y(C—H) 
C=CH;) 
(CNH) , 


ma #7. out-of plane bending: p, 


focoa w, weak 


rte 0 the wavenumber indicates the 


ee eae 
i... at a Se 


IR 


a 
v, 3075(5.3) 
v, 3035(5.0) 
v, 1600(5.0) 
v4 1475(8.0) 
vs 1276(8.8) 


v, 1095(7.9) 
v, 741(8.0) 
vg  686(7.9) 


vg  548(4.9) 
Yio 330(5.6) 
v1, 298(5.8) 


a2 
Vi2 898(4.7) 


v13 498(6.4) 
Viq 168(5.9) 
by 
Vis 862(8.8) 
Vie 765(8.0) 
Vi7 465(4.1) 
Vig 287(6.2) 
Vig 260(5.4) 
Yoo 155(5.9) 


by 
Y21 3022(7.2) 
V22_ 1530(4.2) 


25 1380(6.4) 
Voq 131 8(6.4) ; 


Vos 1 232(6.6) 


tae 1167658) 


Raman Vibration 
number 
2,6-DMA 
-- 2 
-- 20, 
1602(s)p 8, 
1477(w) 19, 
1277(w) 14 
1097(w) 2 
-- 12 
682(s) l 
552(s) 6, 
ae 18, 
292(m) 15 
-- 17, 
500(w) 16, 
170(s)\dp —-:10, 
= 11 
- 4 
287(s\dp_ 10, 
262(s)dp “et 
= 17, 


Sh, ae hes vo 


( bea <a a a 


» | 


en 


ee 
~1) of 2,5-DMA and 2,6-DMA 


Assignment 


v(C—H) 
v(C—H) 
v(C—H) 
v(C—C) 
v(C—C) 
Kekule vibration 


6(C—H) 
d(C—C—C) 
v(C—C) 

Ring breathing 
d(C—C—£) 
6(C—CH;) 
6(C—NH,;) 


(C—H) 
oe 2 


iwe 


os HS age 
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Table 2—Combination and Overtone Frequencies (incm~ ') 
2,5-DMA 
3215 = 3048 + 170 = 3218 
2920 = 1460 x 2 = 2920 
2735 = 1510 + 1220 = 2730 
2620 = 1310 x 2 = 2620 
1800 = 1310 + 490 = 1800 
1720 = 1510 + 215+ 1725 
1702 = 1150 + 550 = 1700 
1423 = 1220 + 200 = 1420 
1296 = 762 + 535 = 1297 


2,6-DMA 
2918 = 1530+ 1380=2910 A, 
2856 = 1475 + 1380 = 2855 B, 
2738 = 1600 + 1145 = 2745 B, 
2632 = 1318 x 2 = 2636 A, 
1888 = 1600 + 287= 1887 B, 
1830 = 1145 + 686 = 1831 _B, 
1770 = 1600+ 168 =1768 4, 
1442 = 1276+ 168 = 1444 A, 
1012 =686+330=1016 4 


863 = 490 + 370 = 860 995 = 498 x 2 = 996 A, 
585 = 370 + 215 = 585 600 = 298 x 2 = 596 A, 
390 = 170 + 215 = 385 575 = 287 x 2= 574 A, 


>~A A RRR RRARR RAD 


325 = 150 + 170 = 320 558 = 298 + 260 = 558 B 


530 = 240 + 287 = 527 B, 


Table 3—Internal Modes of Vibrations of the CH; Group 


2,5-DMA 2,6-DMA Assignment 
2985 2970 vs(CH 3) 
2945 2942 v,s(CH) 
2895 2892 v, (CH;) 
2832 2832 26,,(CH;) 
2710 2710 26, (CH;) 
1370 1375 6,,(CH3) 
1303 1312 6, (CH;) 
1012 1002 p, (CH3) rocking 


Table 4—-Internal Modes of Vibrations of the NH, Group 


2,5-DMA 2,6-DMA Assignment 
3440 3460 v4; (NH) 
3355 3377 v, (NH;) 
1620 1616 v, (NH,) scissoring 
1040 1035 6,; (NH) rocking 
745 740 6,; (NH) wagging 
610 600 y, (NH,) twisting 


3.2 C—C Stretching Vibrations 
The parent molecule benzene has six modes of C—C 
stretching vibrations 8,, 8,, 19,, 19,, 14 and | having 
frequencies 1585 (¢3,), 1485 (¢;,), 3010 (b2,) and 992 
(a,,)cm_ '. The frequency of first five modes remains 
unchanged on substitution, while the frequency of 
mode | (ring breathing) changes in C,,. symmetry”. The 
frequency of the mode 2, changes by a small 
magnitude on substitution in C, symmetry. In the case 
of 2,5-DMA, the bands v4 (R 1638), vs (IR 1578, R 
1588), v, (IR 1510, R 1500), v> (IR 1460, R 1460), vio (IR 
1220, R 1224) and v, (IR 565, R 575) have been assigned 
= as C—C stretching vibrations corresponding to the 
modes 8), 8,, 19», 194, 14 and 1 respectively. In the case 
6-DMA, the bands V3, 22. V4; Y23- Ys and Vs have 
been assigned as C—C stretching vibrations 
corresponding to the modes 8), 84, 19), 194, 14 and 1 
respectively. The frequency of mode | has been 
signed at vg (IR 686, R 682) as ring breathing 


j 
2 
> 
7 
} 
7 


vibration. These assignments are in good agreement 
with the assignments given by Stojiljkovic and 
Whiffen'® and Venkateshwarlu and Radhakrishna?!. 


3.3 C—-H, C—X In-plane Bending Vibrations (4) 

The in-plane bending vibrations of trisubstituted 
benzenes are derived from the C—H in-plane bending 
vibrations of benzene. The six in-plane bending 
vibrations in Wilson’s notations are 3,9,,9,, 15,18, and 
18,. In trisubstituted benzene derivatives, three modes 
of vibrations remain almost unchanged and are called 
C—H in-plane bending vibrations, while three modes 
change considerably in frequency and are called C—X 
in-plane bending vibrations. The modes 3, 9, and 9, are 
regarded as C—H in-plane bending vibrations in 
vicinal trisubstituted benzenes (C,,.), while the modes 
1S, 18, and 18, become C—X in-plane bending 
vibrations. For 2,6-DMA, the bands v6, v2.5 and V6 
have been assigned as 9,, 3 and 9, modes of vibration, 
while the bands v,9 (IR 330) and v3 (IR 240, R 242) 
have been assigned as C—-CH, in-plane bending 
vibration and v,, (IR 287, R 292) have been assigned as 
C—NH, in-plane bending vibration. For 2,5-DMA 
(C,), the bands V9, v,, and v, 3 corresponding to modes 
3, 1S and 18, respectively have been assigned as C—H 
in-plane bending vibrations while v,9 (IR 290, R 290), 
Vo (IR 215, R 215) and v,, (IR 185) have been assigned 
as C—CH,;, C—CH, and C—NH,j in-plane bending 
vibrations respectively. These assignments are in good 
accordance with those of Polyler'? and the frequency 
range agrees well with that given by Varsanyi’. 


3.4 C—C—C In-plane Bending Vibrations 

The normal modes 6,, 6, and 12 are regarded as the 
C—C—C in-plane bending vibrations. For 2,5-DMA, 
the bands v, 5. ¥;7 and Vv, have been identified as C— 
C—C in-plane bending vibrations while in the case of 
2,6-DMA, the bands of v,7 (IR 741), Vo (IR 548, R 552) 
and V5o (IR 455) have been assigned as C—C—C in- 
plane bending vibrations. The above assignments are 
within the frequency range given by Varsanyi> for 
similar compounds and find support from the works of 
Green et al.’ and others®~ '?. 


3.5 Out-of-plane Bending Vibrations (7) 

The six »(C-——H) modes of vibration of benzene are 
5, 174. 17», 10,, 10, and 11 with frequencies 985 cm~' 
(b,), 970 cm — : (@>,)s 898 cm : (4) and 671 cm | A (a>,) 
respectively. Out of these, the a3, fundamental (mode 
11) has been identified as the out-of-plane motions of 
C—H group in phase with respect to each other and is 
highly characteristic! ’. In the case of 2,5-DMA (C,), the 
»(C—H) vibrations will correspond to modes 5, 17, 
and 11 and the three remaining will be »(C—X) 
corresponding to 10,, 10, and 17,. Similarly, in the 


165 


INDIAN J PURE & APPL PHYS, VOL 


compound 2,6-DMA (C;,), the three y(C—H) 
correspond to modes 17,, 17, and 11 and the three 
y(C—X) to the modes 5, 10, and 10. The frequencies 
V>> (IR 802, R 805) anc v1.6 (IR 862) have been assigned 
as the »(C—H) vibrations of mode 11 in 2,5-DMA and 
2,6-DMA respectively. The rest y(C—H) frequencies 
have been assigned at V2, (IR 370, R 375) and V29 (IR 
200) for 2,5-DMA, ane v,> (IR 898) and V9 (IR 155) for 
2,6-DMA respectively. For 2,5-DMA (C,); the 
frequencies V7 (IR 230), ¥29 (IR 170) and v39 (IR 150) 
have been assigned as y(C—CH3;), y(C—CH3;) and 
y(C—NH,) respectively. Similarly, for 2,6-DMA (C;,), 
the frequencies v,4 (IR 168, R 170), v,g (IR 287) and Vy 
(IR 260, R 262) have been assigned as y(C—CH;3), (C— 
CH,) and (C—NH.,) respectively. These assignments 
agree well with those given in the works of Green et al.’ 
and others®~'°. As regards the degenerate pair 16, 
Shimanouchi et al.'* calculated the frequencies of 
asymmetric trialkylbenzene by normal coordinate 
analysis and obtained the values 588 and 441 cm! 
respectively. Bentley and Wolforth'® assigned the 
component 16, between 530 and 590 cm“! and 16, in 
the interval 430-440 cm ' in the case of asymmetric 
trialkyl benzene. In view of these data, the assignments 
of v4 (IR 550) and v,; (IR 535, R 528) in the case of 2,5- 
DMA and v,; (IR 498, R 500) and v,, (IR 465) in the 
case of 2,6-DMA to the modes 16, and 16, respectively, 
may be taken to be justified. The calculated values for 
these modes 16, and 16, for vicinal trimethylbenzene 
(Cz,) are 482 and 487 cm™! as reported by 
Bogomolov!®. Hence our assignments for these modes 
are thus in harmony with the results of earlier workers 
on similar compounds. The normal mode 4 has been 
identified at v,, (IR 730, R 733) for 2,5-DMA while at 
Y16 (IR 765) for 2,6-DMA. Whiffen'’ assigned mode 4 
at 757 cm™' in the case of aniline. Hence our 


assignment for the same mode 4 is in agreement with 
the above data. 


3.6 Internal Vibrations of NH, Group 

Two bands appear in the region 3300-3500 cm~! in 
the spectra of primary aromatic amines associated with 
N—Hstretching vibrations!®:'°. The first band occurs 
near 3500 cm™' due to asymmetric stretching mode 
_ while the second appears near 3400 em ' due to N—H 
symmetric Stretching mode. The asymmetric and 
symmetric stretching vibrations in 2,5-DMA have been 
identified at 3440 and 3355 cm~! respectively. The 
same mode has been identified at 3460 and 3377 em-! 
respectively in 2,6-DMA. The scissoring vibration in 
aniline has been identified at 1618 cm! by Evans?° 
This result confirms the assignments of the frequencies 
1820 and 1616cm~' to the NH 2 SCissoring vibration in 
the molecules 2,5- and 2,6-DMA respectively. The 
bands identified at 1040 and 1035 cm~" in the 
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molecules 2,5- and 2,6-DMA have been assigned pe 
NH, rocking vibration. Cleaves and Dakin? 

identified NH, wagging mode at 783 cm in 
methylamine. In the present case, the frequencies 745 
and 740 cm~! have been assigned as NH, wagging 
mode for the molecules 2,5- and 2,6-DMaA respectively. 
The out-of-plane symmetric NH, stretching vibrations 
have been identified at 510 and 600 cm * for 2,5- and 


2,6-DMA respectively. 


3.7 Internal Vibrations of CH, Group 

In the methyl group there are three C—H stretching 
vibrations, one being symmetric and the other two 
asymmetric. The frequencies of asymmetric vibrations 
are higher than those of symmetric ones. Green et al.’ 
assigned the asymmetric C—H stretching vibrations 
of methyl group in the region 2900-2980 cm~‘' and 
symmetric C—H stretching vibrations in the region 
2820-2940 cm '. In the present case, we have assigned 
frequencies 2985, 2945 cm‘ (for 2,5-DMA) and 2990, 
2942 cm‘ (for 2,6-DMA) as C—H asymmetric 
stretching frequencies. The C—H symmetric stretch- 
ing frequencies have been identified at 2895 and 2892 
cm‘ in 2,5-DMA and 2,6-DMA respectively. In the 
case of 2,5-DMA, the frequencies 2832 and 1370 cm! 
have been assigned as C—H asymmeiric in-plane 
bending vibrations while frequencies 2710 and 1303 
cm * have been assigned as C—H symmetric in-plane 
bending vibrations. Similarly, in the case of 2,6-DMA, 
the asymmetric in-plane bending vibrations have been 
identified at 2832 and 1375 cm ' and the symmetric in- 
plane bending vibrations have been identified at 1312 
and 2710 cm '. Thr rocking vibration of CH3, i.e. p, 
(CH;) has been identified at 1012 and 1002 cm~' in the 
molecules 2,5- and 2,6-DMA respectively. These 
assignments find support from the work of Singh and 
Prasad” and are within the frequency intervals given by 
Varsanyi’. 
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An ultrahigh capacitance ratio silicon hyperabrupt varactor, fabricated by a simple ee ae ae 
Poisson’s equation is solved to derive the expressions for the capacitance-voltage characteristics of t e structure. capac 
ratio, which is defined as the ratio of the diode capacitance at 3.0 V to that at 25.0 V, as high as 9 is obtained. Measurements 
made with the device are compared with corresponding theoretical calculations and an excellent agreement is found. 


1 Introduction 

A varactor diode has found extensive application in 
switching or modulation of a microwave signal, 
harmonic generation and parametric amplification 
and as a tuning element! ”” by virtue of its highly 
nonlinear capacitance-voltage characteristics. It is well 
known that a hyperabrupt P-N junction exhibits a 
large capacitance variation® '* and thus provides a 
higher capacitance ratio for the given two bias voltages 
compared to the conventional abrupt or linearly 
graded structures. For television tuners where varactor 
finds a wide application, the capacitance ratio R is 
defined’’~ '° as the ratio of the capacitances at 3.0 V 
and 25.0 V, ie. R=C;/C,5. In recent years, many 
sophisticated techniques like ion implantation!®!7, 
molecular beam and multi-epitaxy'*~?° or a suitable 
combination of these techniques?! have been used to 
fabricate such high capacitance-ratio varactors. But we 
believe, so far no one has reported a varactor having a 
capacitance ratio more than 6, fabricated by some 
simple technique. Recently, in an attempt in this 
direction, we have reported’? a simple planar process 
wherein ion-implantation in combination with the 
conventional thermal-diffusion techniques was 
employed to fabricate such a device. 

In this paper, we report a silicon varactor having 
very high capacitance ratio of the order of 9, which is 
higher than any Previously published measurement. It 
has been fabricated by a simple chemical vapour 
deposition (CVD) technique in combination with the 
conventional thermal diffusion. The solution of 
Poisson’s equation has been presented here for the 
structure, to describe its C-V characteristics. The 
fabrication procedure, C-V characteristics and 
impurity Profile of the device are presented and 
compared with the theoretical results. 


2 Theory 
A simple expression for capacitance law of a 
hyperabrupt structure can be obtained when the 
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carrier distribution in the active region of its junction is 
only of one type, i.e. Gaussian”? or exponential!?. In 
the literature numerical solution of the Poisson’s 
equation has been reported!°'!* to describe the C-V 
characteristics of a hyperabrupt junction which 
contains more than one type of distribution in its active 
side. No simple C-V expressions are available for such 
structures. Here we have derived the C-V relations fora 
hyperabrupt structure in which active region consists 
of a Gaussian type retrograded region formed on a 
uniformly doped epi-layer as shown in Fig. 1. The other 
side of the junction is P-region which is doped heavily 
to avoid depletion of carriers in it. 

To describe electrical characteristics of the structure. 
we have to solve one-dimensional Poisson’s equation 
for two regions, viz. (1) for X < X, and (ii) for X > Xs; 
with the boundary conditions: 

V,(0)=0 
V(X) = V2(X4) 
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Fig. 1—Typical impurity profile of ultrahi capacitance rati 
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dV,/dX\x=x, =dV2/dX|y~ x, 
V,(W)=(V,+ V;) 
dV,/dX |y-y =0 


In the above expressions, V is the electrostatic 
potential and subscripts 1 and 2 refer to first and 
second region respectively, W is the depletion layer 
width, V, is the applied reverse bias. V; is the built-in 
potential and X, is the width of the retrograded region. 
The carrier distribution in the active region is given by 


N(X) = N, exp(— X?/a?) i 

oo for XS X, “ 

fim 5 ae, & 
y= Ne ae oy 


Here ‘N ,’ is the concentration at the junction X = 0, 
a the diffusion constant and Nz, the background 
concentration. Using the boundary conditions (1) and 
the distribution given by Eqs (2) and (3), the solution of 
Poisson’s equation gives 


I 


= |  _______— sak 
= Saeiag | 2KeolVa + V;) *) 
qN,a? 
x 
fy for X< X, 
2K é.(V,+ V;) N, 
W = ed at (Ft 1) 433 | 
J qNp Neg ot tS) 
for X > X, 


where K is the relative dielectric constant of silicon and 
Eo is the permittivity of free space. 

The junction capacitance is given by 
C=Ke,A/W ... (6) 
where A is the active area of the junction. Using Eqs (4) 
and (5) in Eq. (6), the C-V characteristics of the 
structure can be described as 


Ké A 


C= l 
‘ lm ara | 


qN ,a? At 
for the applied bias for which WS X, and 


KéQ A 
2K &)(V, + a ae at 


| qNp 


for the bias for which W> X,. | 7 

It is assumed here that the thickness of the epi-layer ts 
large enough so that the depletion layer edge does not 
penetrate into the out-diffusion tail of the epi-substrate 
for the maximum applied voltage. If N, and a are such 


that the depletion layer confines within X , for V,=3V 


C= 


ae 


and enters into the epi-region for higher bias voltage, 
the capacitance ratio R is given by the ratio of Eq. (7) 
and Eq. (8) for V, = 3.0 and 25.0 V respectively. It is 
obvious from the ratio of Eq. (7) and Eq. (8) that with 
suitable choice of N,, Ng and a, a large value of R can 
be achieved, and for this a large value of N, and small 
value of Nx are preferred but their values are restricted 
by the requirements of high breakdown voltage** and 
low series resistance'?''* of the device respectively. 


3 Diode Structure Fabrication 

The diode structure is shown in the inset of Fig. 1. 
The retrograded region is formed in n/n* epitaxial 
wafer by chemical vapour deposition (CVD) technique. 
The silicon wafer containing diffusion windows in 
silicon dioxide is coated with phosphorus-doped silox 
layer, grown by the reaction of oxygen and silane 
containing phosphine at 450°C. It is then heated to a 
high temperature in nitrogen ambient to obtain the 
desired depth X,. The junction is formed by diffusing 
boron into it using BN source. Subsequently silicon 
oxide is grown by thermal method and the contact 
windows are opened. The wafer is then aluminized and 
the devices are encapsulated in TO-18 headers. 


4 Results and Discussion 

The C-V characteristics and the impurity profile of 
one of the diodes fabricated as above are shown in Fig. 
2. C-V measurements were carried out by using 
Boonton Capacitance Meter Model 72 AD, and the 
profile was computer calculated using the measured C- 
V data. The package capacitance which is the 
capacitance of an open diode in the package has been 
subtracted in the above measurements. From Fig. 2 it is 
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Fig. 2—C-V characteristics and impurity profile (N versus W) of 
ultra high capacitance ratio silicon varactor. 
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seen that the capacitance ratio is 9 and the width of the 
retrograded region is 1.1 um. 
Calculation of retrograded impurity profile 
described by Eq. (2) gives N; = 4.8 x 10!° atoms/em* 
and a=6.2x10~° cm. These values give the best fit 
with the measured profile, as shown by triangular 
points marked in the Fig. 2. Eq. (4) gives V, = 12.6 V for 
W=1.1 pm. C-V characteristics are calculated by Eq. 
(7) for V,=0 V to V, = 12.5 V and by Eq. (8) for V, 
> 12.5 V. Built-in potential V; = 0.8 V is taken in the 
calculations as measured by sie of 1/C? versus bias. It 
can be seen from Fig. 2 that an excellent agreement 
exists between experimental results and theoretical 
calculations. Eqs (7) and (8) give the value of R= 8.8 
_ which agrees well with the average value 9 given by Fig. 
_ 2. Thus the derived expressions (7) and (8) describe the 
_ C-V characteristics of the structure fairly accurately. 
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A simple phase shifter developed for the reference signal is capable of handling the entire range of field modulation 
frequency used in magnetic resonance spectrometers. A discussion of the normally adopted technique for the reduction of the 
unavoidable switching spikes which arise when JFETs are used as analog switches in the phase-detector, is also given. 


1 Introduction 

Phase-sensitive detectors. yield the ultimate 
information about the height and width of the spectra. 
The signal-to-noise (S/N) ratio. of the detector output is 
poor for certain samples. So the design of the phase 
sensitive detector is rather critical in the entire 
instrumentation of the magnetic resonance spectro- 
meter. Description of such phase-sensitive detectors 
using different circuit techniques have appeared in the 
literature’ *. Improvement in the S/N ratio is 
achieved with the help of a gradual phase shifter for the 
reference channel which does the job of maximization 
of the phase detected output. The full range of the 
phase shifter needs to be 0° to 360° with an additional 
constraint of the minimum variation of the magnitude 
of the phase shifter output for any arbitrary change of 
phase. These two conditions cannot be met at high 
modulation frequencies without the use of active 
devices. Phase shifter circuits with active devices 
consist of a paraphase inverter and a series R-C 
combination between the two outputs of the paraphase 
amplifiers*. This method suffers from the drawback 
that small phase shifts in the neighbourhood of 0° and 
180° are not obtained in actual practice due to an 
additional phase shift caused by the coupling capacitor 
and the load resistors of the paraphase amplifier. So, 
incorporation of two such phase shifting circuits does 
not lend itself to any arbitrary phase shifting of the 
reference channel between 0° to 360° which may be 
required for the maximization of the detected signal. 


Another technique of phase shifting using an 


operational amplifier with R-C network has been in 
use for signal processing applications. Here also, the 


same back of the discontinuity of the gradual 
Bee at neighbourhood of 


audio frequency range. So, this is also not adaptable for 
the present purpose.’ 

A digital method of phase shifting’ has been used in 
some commercial instruments which is somewhat 
more complex and requires a locked second harmonic 
signal with respect to the sinusoidal modulation 
frequency. This method derives a double frequency 
square wave through zero-crossing detector and 
utilizes this to achieve 2 phase square waves at 
fundamental frequency. These signals have phase shifts 
0, 90, 180 and 270° with respect to the field modulating 
sinusoidal wave. A digital method for quadrant 
selection is used and the same principle of R-C phase 
shifting for analog signals is used with the help of 
interface circuits where the analog voltage across the 
capacitor is fed to a comparator for bringing back the 
phase shifter output in the form of logic signal. This 
logic signal is used as clock for triggering a flip-flop 
wherefrom complementary pulses are available for the 
phase detection. These pulses are phase shifted with 
respect to the square wave obtained from the reference 
channel through a zero-crossing detector. The method 
features excellent performance but involves a large 
number of components which are used in the 
cumbersome combination of analog and digital 
circuits. 

The methodology presented in this paper claims the 
novelty of an analog phase shifter which is capable of 
providing a continuously variable phase output, in 
relation to the reference sinusoidal signal and yet is 
simple in its configuration. 


2 Design of the Pre-amplifier 

The pre-amplifier circuit for a field-modulation 
frequency of 100 kHz is presented in Fig. 1. The front- 
end amplifier consists of an R-C coupled amplifier with 
JFET as the active element to ensure a high input 
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impedance at the frequency of operation. The voltage 
gain is further achieved with a single stage of bipolar 
transistorised R-C coupled amplifier. The frequency 
selectivity is obtained with the help of an acceptor 
circuit using passive components and is being driven 
from a buffer amplifier fed from the output of the 
bipolar stage. This buffer amplifier is a power amplifier 
where a superalpha device, comprising an N channel 
JFET and a PNP power transistor, is used as the active 
element. Normal power amplifiers employing bipolar 
junction transistors as impedance converters would 
oscillate at this frequency due to improper isolation 
between the input and the output. The JFET in the 
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superalpha combination acts like a source follower and 
provides a good isolation between output and input 
while the PNP transistor, connected in common 
emitter mode, is directly coupled to the JFET stage 
operates without the slightest tendency of oscillation. 
This combination provides a linear output with a very 
low impedance so as to deliver the current through a 
series resonant circuit. The series resonant circuit 
presents a constant resistive impedance of 120 ohms 
under resonance. The attenuator is connected in 
cascade with the resonant circuit such that the 
impedance offered to the resonant circuit does not 
change with different output voltage level. The 
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attenuator output drives an amplifier being made out 
of a fast operational amplifier under simple inverting 
configuration. The amplifier power level is then 
boosted with the help of a superalpha combination as 
described earlier. The transformer is used for the 
purpose of obtaining electrically-isolated push-pull 
signals which is required for the phase detection. This 
transformer is being driven by the output of the power 
amplifier made of superalpha pair. A simple reactive 
compensating network is added in parallel with the 
transformer primary so as to present an overall nearly 
resistive load to the amplifier. 


3 Phase Sensitive Detector 

The basic purpose of phase-sensitive or lock-in 
detection is the retrieval of the signal with a high S/N 
ratio. In order to maximize the signal, the phase of the 
carrier has to be adjusted keeping its magnitude fairly 
constant. The present scheme adopts a novel circuit for 
continuous alteration of the phase at any frequency 
with R-C networks and active devices, and does not 
require any special components. The circuit 
configuration as presented in Fig. 2 is uniquely 
versatile in respect of its adaptibility over the entire 
range of field modulation normally employed in 
magnetic resonance spectroscopy. 

The basic phase shifter consists of an IC differential 
amplifier (CA 3028A) whose collectors are connected 
through a series R-C combination so that the phase 
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shift at the junction point of the resistance and 
capacitance is exactly 90° with respect to the input to 
the differential amplifier. This is again fed to another 
differential amplifier (CA 3028A) with unity gain so as 
to obtain 90 and 270° phase shifts with respect to the 
original signal, i.e. the carrier. So the four outputs of the 
differential amplifiers will have same amplitude with 0, 
90, 180 and 270° phase shifts with respect to the carrier. 
A potentiometer connected between 0 and 90° will give 
an Output voltage with variable phase. Similarly, phase 
variation in other quadrants is also achieved. Thus the 
detected signal is maximized by changing the phase of 
the reference signal. 

The magnitude variation of the phase shifted voltage 
is not zero as One expects with this type of circuit. 
However, from the following simple analysis the 
estimate of magnitude variation can be obtained. 

Let the two signal sources connected to the two 
extreme points of the potentiometer be E L0° and 
E £90° with a common ground reference. The total 
resistance of the potentiometer is R. For a wiper setting 
of a R and (1 — a)R to the left and right of the wiper, the 
voltage at the wiper with respect to common reference 
is given as 
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The value of |Eo| |E| is observed to vary from unity to 
0.707 for the entire range of values of a between Oand I. 

This variation of voltage does not impair the 
detector operation due to the fact that switching 
signals have been derived from the output using Zero- 
crossing detector (710 comparator). This output is used 
for driving widely used JFET switches as described in 
the following paragraph. 

The output from the phase shifter is fed to a super- 
alpha pair for obtaining the low impedance required 
for the following stage of amplification using video 
amplifier type 733. The differential outputs from this 
amplifier are fed to the inputs of the analog comparator 
710. Due to the differential signals, the jitter in the 
Output of the comparator, used without positive 
feedback, is totally absent. This comparator output is 
used to trigger a monostable multivibrator 74121 
whose delay time is adjusted with external R-C 
components so that the square wave at the output is 
exactly symmetrical and is independent of the input 
conditions. The complementary outputs are level 
shifted using PNP transistors 2N 2905 to derive the 
conventionally used JFET switches. The usual 
technique® of the suppression of the switching spikes 
occurring at the output of the phase detector is also 


Fig. 4—(a) Typical ESR spectrum of Mn?* ; 
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incorporated in this design. An attempt has been made 
to explain this spike reduction mechanism in the 
following paragraph. 

The change of the depletion layer charge during a 
switching operation on a JFET is solely responsible for 
the spike at the channel. If the gating signal to the JFET 
is coupled through another JFET, the problem can be 
overcome. successfully. Then the limitation on the 
frequency of switching is only due to the transit time 
phenomenon and the dominance of drain-source 
capacitive reactance. 

The basic scheme of the JFET switch with a JFET- 
connected diode is shown in Fig. 3(a). The resistance R 
merely provides a gate leak. so that in the event of the 
diode-connected JFET being under cut-off condition. 
the gate may not pick up any stray noise. Fig. 3(b) 
shows the depletion layers of both the JFETs when the 
switch is Off. The diode-connected JFET becomes 
forward biased.so that the depletion layer shrinks toa 
minimum width and this entire depletion layer charge 
is being poured on to the depletion layer cf the 
switching JFET. causing a pinch-off in the channel. 
Similarly, while the switching JFET is On, its depletion 
layer width is minimum and the entire charge in its 
depletion layer is now accumulated into the depletion 
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layer of the diode connected JFET which is now under 
cut-off condition. This is shown in Fig. 3(c). Therefore, 
the pair of JFETs looks like a combination of a charge 
source and a charge sink, and the charge moves from 
one JFET to the other, depending on the gating signal. 
Consequently, there is no charge transfer external to 
the switch and its driver combination. As a result, this 
_ JFET pair may be used as a switching device up to a 
_ few megahertz of frequency as the input capacitance 
and Miller capacitance are both very small for JFETs. 
_ Circuits following the detector are a dc amplifier 
(OPO7), second-order low-pass active filter (LH 0042) 
and a recorder pre-amplifier for signal recording. 


4 Performance 

_ The delay introduced by the monostable multivib- 
ator (one-shot) (74121) depends on the external 
Seagle whose aging and instability are the only 
mM 2 points against this design. However, the 
y stability of the reference modulation signal 


1 (minimum is | part in 10°). The percentage 
ior of the delay time set to produce symmetrical 
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insignificant for the overall performance of the phase 
detector output. 

The performance of the phase detector and the phase 
shifter can best be checked with regard to its stability, 
drift and resolution of the signals by taking an ESR 
spectrum of Mn?” in calcite with Mn content ~ 27 
ppm (Fig. 4a) as well as colour centres in the yellow 
variety of CaF, (Fig. 4b). 
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An ultrasonic viscometer for the measurement of dynamic shear viscosity of liquids at 10 MHz frequency over a wide range 
of temperature has been described. Shear wave reflectance technique at normal incidence has been adopted. A BT-cut quartz 
crystal along with a fused quartz delay line has been used for generating shear waves. The viscosity measurements performed on 
glycerine reveal that the results are in good agreement with known accepted values. The viscometer has been found to be 
capable of measuring viscosity with comparable accuracy under temperature-controlled conditions. The system is reliable and 
offers potential applications for studying viscoelastic behaviour of liquids in the frequency range (3 to 50 MHz) by 


incorporating some minor modifications. 


1 Introduction 

Ultrasonic techniques'’? are known for measure- 
ment of dynamic viscosity of liquids. This parameter is 
required to be known in industrial applications such as 
lubrication, where the liquids are subject to alternating 
shear rates***. The use of torsional piezoelectric 
crystals*’® for measurement of viscosity is limited up to 
a few poises only at low ultrasonic frequencies in the 
kHz range. The few commercially available 
viscometers’ based on this technique are not suitable 
for finding the dynamic viscosity of liquids in the high 
ultrasonic frequency range. 

In studying the viscous behaviour of liquids in high 
frequency range 3-250 MHz, though the ultrasonic 
shear reflectance technique®~!° has been successfully 
employed, no compact instrument is presently 
available. The need for the development of a suitable 
viscometer is felt and it can be expected to have a wide 
scope in various industrial applications. 

An attempt is, therefore, made to develop an 
ultrasonic viscometer at the National Physical 
Laboratory, New Delhi, using the normal incidence- 
reflectance method. The instrument has been 
successfully used for the measurement of dynamic 
viscosity of glycerol in the temperature range 4-33°C. 
In this paper, the ultrasonic viscometer developed by 
the authors is described along with its working 
Principle, and the procedure for measurement of the 
dynamic viscosity of liquids outlined. 


2 Theoretical Considerations 
The shear viscosity of a liquid under investigation is 


determined experimentally by the measuremeat of the 
reflection coefficient of an ultrasonic shear wave 


incident at the interface of the liquid and a fused quartz 
ultrasonic delay line of known shear impedance. A BT- 
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cut quartz transducer fixed on one end of the fused 
quartz rod generates pulses of ultrasonic shear wave, 
which travel along the delay line and is reflected at the 
interface. The complex coefficient of reflection K is 
given by the expression®: 


K =(Z,, — Z)/(Z 54 + Z,) = kexp(— j6) ey 


where k is the absolute value of the reflection coefficient 
and @ the phase shift. Zyq=Rygt+ jX sq and Z,=R, 
+ jX;, are the shear mechanical impedances of the 
fused quartz and the liquid respectively. 

Using Eq. (1), R, is given by the relation® 1 = 


R, = Ry,(1 —k)(1 +k) at 


Where Ry, is known® to be equal to 8.3x 105 
mechanical ohm (dyne-s/cm). Now the magnitude of 
the reflection coefficient (k) can be evaluated by the 
expression!!; 


k =(Am+n/Am)iq/(Am +n/Am)air saan 


where A,, and A,,,, are the amplitudes of the mth and 
m+ nth echoes respectively. The dynamic viscosity 9 is 
the viscous component of the shear modulus and is 
given by the known relation?**:°: 


n=R?/(pxf) ...(4) 


where p is the density of the liquid and f the frequency 
of shear waves. 


3 Constructional Details 
The viscometer consists of an ultrasonic measuring 
cell which contains a quartz crystal, quartz delay line, — 
and liquid container. Its electronic system comprises a _ 


an oscilloscope. The block diagram of 
System is shown in Fig. 1. 
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visualization On the screen. The schematic pattern of 
oscillogram of the received signal is shown in Fig. 3. It 
is a series of peaks of decreasing amplitude which 
represent the echo intensity after successive reflections 
from the interface. The magnitude of the reflection 
coefficient is found out by substituting the measured 
amplitudes of the echoes in Eq. (3). 

The Ultra-Cryostat type MK-70 has been used for 
measurement of viscosity of liquid under investigation 
at various temperatures ranging from 33°C to 4°C. 
However, it is capable of maintaining the temperature 
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QUARTZ USED from + 60°C to — 30°C. The measuring cell containing 
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temperature bath. 


4 Measurements and Discussion 
To check the performance of the viscometer 
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~ — ; ° 4° 
Table 1—Measured Values of Shear Viscosity of Glycerine Over the Range of Temperature 33°C to 4°C 


Temp Reflection R, Value of np Frequency Ref. No. - S : 
te coefficient § 10° mech.Q P-gem ~* MHz g 
k 
Present Literature 
set-up value 
3.0 0.97 12.60 5.08 5.01 14.7 (13) 1.252 4.05 
- 0.96 17.60 9.85 11.32 14.7 (13) as Bis 
0 0.95 21.18 14.41 —_— — _ : n 
: 0 0.93 29.08 26.89 28.10 0.5 (12) 1.268 21.20 
7.0 0.92 36.07 41.41 42.52 0.5 (12) «1.270 ~—- 32,60 
4.0 0.89 48.30 14.27 72.50 0.5 (12) (1.273. ~—«$8.34 


standard liquid for calibration purposes, it was chosen 
merely because of its being a well-studied liquid and 
accepted values of viscosity over wide ranges of 
frequency and temperature are available in litera- 
ture’*~**. The reflection coefficient of shear waves at 
fused quartz/glycerine interface was measured at 10 
MHz in the temperature range 33-4°C. This, in turn, 
was used for determining the real part of shear 
impedance and hence its viscosity using Eqs (2) and (4). 
Measured values of the product of viscosity and 
density of glycerine are summarized in Table 1 along 
with the results obtained by earlier workers for 
comparison. 

The variations of np values with temperature as 
observed by the different workers are compared in Fig. 
4. It is seen from the Fig. 4 that the np value for 
glycerine increases from 5.08 to 74.27 P g/cc as the 
temperature is decreased from 33°C to 4°C. The results 
obtained with the set-up developed are found to be in 
agreement with the earlier reparted results within 4 Te 

Assuming the density values of glycerine from the 
Critical Table'*, the values sof viscosity were 
determined from np values. It is seen from Table 1 that 
the viscosity of glycerine increases from 4.05 to 58.34 Pp 
as the temperature decreases from 33-4°C 


5 Salient Features 

The viscometer described here is reliable and simple 
to Operate over a wide range of temperature. The 
technique can be employed with advantage for taking 
quick sample measurements of viscosity of liquids up 
to three orders of magnitude with an accuracy of 3-5 %/ 
under temperature-controlled conditions. The system 
can also be used for studying the elastic behaviour such 
as real and imaginary components of Shear modulus of 
liquids. Liquid samples of only a few millilitres are 
needed and there is no danger of any depolymerization 
of liquids due to irradiation by ultrasonic waves 
because of the low intensity ultrasound used. 
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6 Concluding Remarks 

The instrument developed at NPL is not only useful 
as a viscometer but also can be used as an impedometer 
for measuring shear modulus of liquids. The system 
will find uses in studying viscoelastic properties of 
polymers and lubricating and machine oils under 
rapidly alternating shear stress. Such conditions are 
normally encountered by lubricating films in high 
speed bearing and gear teeth. 

The accuracy of 3-5 % of the ultrasonic viscometer is 
due to the instrumentation involved in the 
measurement of real component (R,) of the impedance. 
In this system, R; is evaluated from the measurement of 
peak amplitude of the echoes. The accuracy of this 
order is acceptable in routine measurements of 
dynamic shear viscosity of liquids such as lubricants. 
However, the system is under modification for higher 
accuracy by introducing an auxiliary channel for 
measuring the attenuation (in dB) per reflection. The 
instrument can be used for measurement of viscosity of 
liquids at different ultrasonic frequencies by 
incorporating minor modifications in the ultrasonic 


unit and switching over to associated electronic 
circuitry. 
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It is pointed out that the logarithmic potential for short-range 
repulsive interactions in alkali halides used by recent workers does 
not give good accotnt of the cohesive energy. The suitability of this 
potential with and without including van der Waals interactions has 
been examined. The potential parameters and from them the 
cohesive energies are calculated and compared with the éxperimental 
values. It is shown that the logarithmic potential form is inferior to 
the already known and frequently used exponential potential form. 


Recently, Jha and Thakur' have proposed a 
logarithmic potential form for repulsive interaction in 
alkali halide crystals which can be expressed as 
follows: 
w(r) =aln(1 + br~") zeata) 
where r is the interionic distance, and a,b and nare the 
potenital parameters. Jha and Thakur have treated a,b 
and n as free parameters, i.e. their values differ from 
crystal to crystal. However, there exists an error in the 
method of Jha and Thakur? Several similar potential 
forms for certain fixed values of n such as n = 4 or n=9 
have been proposed by previous workers?:3. In the 
present note we examine the suitability of potential 
given by Eq. (1) with and without including the van der 
Waals (vdW) dipole-dipole and dipole-quadrupole 
interactions. When we include the vdW interactions, 
the total short-range potential becomes : 
—& D 
W()=aln(1 + br“) ...Q) 
. of the vdW coefficients C and D are taken from 
the recent calculations based iation. 
ara, on the variational 
For calculating the potential 
normally makes use of the cr 
condition and bulk modulus. Since in the present case” 
there are three unknown parameters, we have to use 
also either the experimental values of cohesive energy5 
or the pressure derivative of bulk modulus®-8 for 


parameters, one 
ystal equilibrium 


tJha and Thakur have taken do(r)/dV = br? (condition (iii) in Ref. 1] 
which is not correct in view of the crystal equilibrium condition. 
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calculating the potential parameters. We have 
calculated the three potential parameters using the 
crystal equilibrium condition, bulk modulus and its 
pressure derivative. This method is similar to that used 
by Michielsen et al.° The following equations are used 
in calculating the parameters 


de(r) fs a 
eee (3) 
(“A”) _ =o highs (4) 
d* g(r) 3 dB 
eis. 2 27KBr( 1 a -.(5) 


where ro is the equilibrium interionic distance, Br is the 
isothermal bulk modulus. dB;/dP represents the first 
pressure derivative of By and g(r) is the total crystal 
lattice energy, written as follows: 


+ W(r) .. (6) 


where a, is the Madelung constant and y(r) is given 
either by Eq. (1) or (2). The potential parameters 
calculated from Eqs (3-5) have been used to estimate 
cohesive energies with the help of Eq. (6). The cohesive 
energies thus calculated are compared with 
experimental values in Table 1, 


Oy C7 
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Table 1—Comparison of Calculated & Experimental Values 
of Cohesive Energies (in kJ/mol) 


Crystal Calc. values Exptl values 
(Ref. 5) 
Eq. (1) Eq. (2) 
LiF 946.5 951.2 1032.2 
Licl 772.4 772.6 850.2 
LiBr 731.5 730.2 812.5 
Lil 654.3 653.9 754.3 
NaF 882.5 881.9 918.4 
NaCl 744.0 743.1 782.8 
NaBr 707.6 708.3 746.8 
Nal 656.6 657.3 698.7 
KF 789.2 788.8 813.0 
KCI 681.2 682.2 7121 
KBr 652.9 654.6 682.8 
KI 614.2 614.8 642.6 
RbF 751.6 750.3 777.4 
RbCl 657.2 659.1 684.5 
RbBr 628.3 * 630.5 657.7 
RbI 593.4 596.8 621.3 
oi 708.6 702.5 740.6 
627.0 636.0 
CsBr 605.2 617.3 


$71.3 


NOTES 


We note from Table | that the calculated values of 

yhesive energies deviate significantly from the 
experimental ones by about 5% to 13%. The 
maximum deviations are in case of lithium halides. The 
inclusion of vdW interactions does not improve the 
agreement between calculated and experimental 
values. On the other hand, it is well known that the 
cohesive energies calculated from the Born-Mayer 
exponential potential’®'' present good agreement 
with experimental values within 1-2%. It may, there- 
fore, be concluded that the logarithmic potential 
is inferior to the well known potentials reviewed by 
_ Tosi’® and used by others at various times. 
___ One of us (RKG) is grateful to the CSIR, New Delhi 
for the award of a junior research fellowship. 
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Curie temperatures for slow cooled and quenched copper ferrite 
samples are used to compute cation distribution using Gilleo’s 
formula. Cation distribution is discussed on the basis of oxygen non- 
stoichiometry and square bond formation of Cu?* in copper ferrite. 
The deviation in theoretical and experimental magnetic moment 
values may be due to randomly distributed Cu** ions in oxygen 
interstices and formation of Fe?* ions on B site. 


Cation distribution in ferrites is usually: found from 
various studies on X-ray diffraction, Méssbauer effect, 
magnetization and Curie temperatures. Recently, 
cation distribution from Curie temperature has been 
reported in Cu-Zn (Ref. 1) and Mg-Mn (Ref. 2) mixed 
ferrites. Curie temperature, which is easy to determine, 
gives the distribution of cations accurately. In this note, 
an attempt has been made to evaluate cation 
distribution in slow cooled and quenched copper 
ferrite samples.-Cation distribution in these samples is 
discussed on the basis of crystal distortion and oxygen 
non-stoichiometry. 

Copper ferrite samples were prepared by the usual 
standard ceramic technique using AnalaR grade oxides 
of CuO and Fe,O3. The details of prepareation, mode 
of quenching and determination of Curie temperature 
were reported elsewhere*. The magnetization of the 
samples was measured on high field loop tracer Hs 869 
(ECIL make) (Ref. 4). 


ee 


The cation distribution is evaluated from Gilleo’s 
formula using the values of Curie temperatures of these 
samples (Table 1). Based on these cation distribution 
values, magnetic moment (/z, in Bohr magneton) is 
computed assuming Neel’s two sub-lattice model. 
From the magnetization studies, the values of ng are 
evaluated using the formula 


Mol. weight x Saturation magnetization 
5585 x Density 


ng= 


It is seen from Table 1 that Curie temperature 
decreases with the increase of quenching temperature. 
This is due to the migration of more Cu?* ions from B 
site to A site. The cation distribution also shows the 
same trend. Rezlescu er al.° have also observed the 
lowering of Curie temperature in Cu-containing 
ferrites. From Table 1, it is also observed that the 
theoretical and experimental values of ng differ and this 
difference is more pronounced at higher quenching 
temperatures. 

The role of oxygen deficiency in changing the cation 
distribution has already been recognized in copper 
ferrite°. Oxygen non-stoichiometry increases both 
with increase in quenching and sintering temperatures. 
The formation of Cu?* to Cu* ions which get trans- 
ferred to A site is favoured by oxygen deficiency. Oba- 
yashi and lida’ have considered the oxygen deficiency 
(0) with quenching temperature and obServed that 
the deviation between the calculated and experimen- 
tal values of ng is more at lower temperatures of 
quenching. They remarked that the deviation in values 
below 900°C may be due to thermally-induced 
distribution of Cu’* ions on A site or the inaccuracy in 
taking the values of 6 (i.e., CuFe,O,-;). The samples 


Table 1—Values of Quenching Temperature, Curie Temperature, Magnetic Moment (ng) and Cation Distribution in Slow 
Cooled and Quenched Copper Ferrite Samples 


Quenching Curie ng (in Bohr magneton) Cation distribution 
temp. temp. calculated from 
Cc °C 
Magneti- Theore- Observed A Site B Site 
zation tical (Ref. 14) 
studies I 
Slow cooled 435 
- % i 2 i 8 = Fe35;Cud to Fe} b9Cud 3, 
600 402 C C 
= i 52 2.92 _ Fe3$¢Cu2 3, Fe} 34Cud 5, 
800 388 1.81 32 ee c Cc 
2 — Fe§$,Cu3$ Fe?5 5 
om io P 2 28 “ef eCuds, 
= 6 08 47 a Fed ieCuas, Fe? 3,Cud is 


studied were prepared at 950°C sintering temperature. 
Oxygen vacancies may not be predominant and 
reduction of Cu** to Cu* may not be possible at this 
sintering temperature®. 

Goodenough and Loeb? have stated that Cu2* ions 
form square bonds on B site in ferrites. The square 
bond formation decreases with increasing quenching 

temperature, resulting in the decrease of tetragonality. 
Finch et al.'° have given the cation distribution at 
which tetragonal to cubic transition takes place in 
copper ferrite as 
7 Cus 34 Fed 66[Cud.66 Fe} 34] 
2 The cation distribution evaluated by Géilleo’s 
_ formula for the sample which is quenched at 900°C is 


Cus. Fedés [Cus 6s Fei3, 
which has also cubic structure found by diffractometer 


+ 


tudies. Bertaut'’ has observed that the cubic state is a 
Stically disordered state. The same feature has 
n observed by Evans et al.'* They concluded that 
he random distribution of Cu** ions in oxygen 
ices in Cu ferrite is more in cubic state than in 
‘State. 

mparison of ng values reveals that when Cu2* 
‘no magnetic moment (i.e. ng=0), the 
of Fe** ions on B site shows satisfactory 
d ng values become very close. Murthy et 

stated that the presence of Fe?* ions 
out even if they could not be detected 
because of their small 


NOTES 


that it gives direct information about cation 
distribution from the knowledge of Curie tempera- 
tures. Work is in progress to generalize this formula for 
ferrites having magnetic ions. 

The author wishes to thank Prof. R N Patil. Head, 
Physics Department, Shivaji University, Kolhapur, for 
encouragement and useful suggestions in this work. 
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The unit cell and space group of heratomin and Hg-derivative of 
geraniol have been determined by X-ray diffraction method. The 
compounds are monoclinic with space groups P2,/m and P2/c 
respectively. 


The work reported in this note forms a part of the 
programme undertaken in this department to 
investigate the crystal structure of some organic 
compounds like C,,;H,40, C,,H,,0;Hg etc. 

Heratomin occurs in Heracleum thomsonni, a herb 
growing in Ladhak area of Jammu and Kashmir state. 
The compound was isolated from the petroleum ether 
(60-80) extract by chromatography over alumina. The 
chemical structure of this compound established by the 
IR, UV, NMR and mass analysis is 6-(3 methyl but 2- 
anyloxy) angelicin’. Crystals of heratomin were grown 
by slow evaporation of the solution of the product in 
acetone-hexane at-room temperature. 

Geraniol is a monoterpene alcohol and its Hg- 
derivative was isolated as an intermediate product and 
its structure was established by the UV, IR, NMR and 
mass analysis methods’. Crystals of Hg-derivative of 
geraniol were grown by slow evaporation of the 
solution of the product in a mixture of benzene and 
methanol (1:1) at room temperature. 

Unit cell dimensions were determined from 
oscillation and Weissenberg photographs taken about 
crystallographic axes with copper radiation filtered 
through a nickel foil. Weissenberg reflection data up to 
Sth layer were recorded for crystals of heratomin and 
up to 4th layer for Hg-derivative of geraniol. The 
reflection data were collected about two crysta- 
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Table 1—Crystal Data of the Compounds Studied 


Heratomin Hg-derivative 
(C,6H,40) of geraniol 
(C,,H2,0;5Hg) 
Molecular weight 270 446 
Measured density 1.32 2.43 
(D,,, g/cm?) 
Calculated density 1.33 2.45 
(D., g/cm”) 
Point group 2/m 2/m 
Space group P2,/m P2/c 
No. of molecules per 8 4 
unit cell (Z) 
a(A) 22.24 17.04 
b(A) 8.92 7.60 
c(A) 14.54 9.33 
V(10- 24cm?) 2691.15 1206.61 
Solvent Mixture of Mixture of CCl, 
CCl, and and tetrabromo- 
hexane ethane 


llographic axes in each crystal. The reflection data were 
indexed by the use of Weissenberg template. The space 
groups for heratomin and geraniol were assigned from 
the systematic absences in the data. The densities of the 
two samples of the compounds were measured by the 
flotation method. The crystal data for the two 
compounds are given in Table 1. Their structure 
determination is in progress. 

We are thankful to Drs S K Banerjee, M K Aggarwal 
and R K Thapa for supplying the crystals and for many 
useful discussions. One of the authors (DR) is grateful 
to the University Grants Commission, New Delhi for 
the award of teacher fellowship tenable at the Jammu 
University. 
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A modified valence force field has been derived for N,N- 
dichloroacetamide from vibrational frequencies using damped least- 
squares method. The nature of the normal vibrations has been 
investigated with the help of potential energy distributions. 
Arguments have been presented to revise the assignment of 
stretching mode of C—C to 1001 cm™!. 


The vibrational assignments of N,N- 
dichloroacetamide (DCA) were made earlier by Devia 

and Carter’ using Urey-Bradley (UB) potential 
function. But the assignment of the band at 1191 cm™! 

| to C—C stretching mode is abnormal and seems to be 
4 ‘in error when compared to the assignment of the 
corresponding band in related systems, such as N- 
chloroacetamide’, acetamide”, N-methylacetamide? 
and N,N-dimethylacetamide* near 955, 874, 894 and 
957 cm”! respectively. Their own calculations, using 
UB model, predicted this mode around 1013 cm ‘'. 
eer itmted this large difference of about 180cm' 
between the calculated and observed frequencies of 
—C) to the omission of significant interaction 


| ¢ twee en shes C—C bond and the methyl group in the — 
JB However, they did take into account all. 


ction constants of significance 


Min cainey 


" wine ae eden: ws 


= 


Cee 
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: -—C bond and the methyl group as they 9 


mode. We have chosen to use a modified valence force 
field, instead of the UB model, for the study of DCA to 
supplement our work of a similar nature on halogen 
substituted amides to understand the effect of halogen 
substitution on the amide characteristic fre- 
quencies®”’. A further advantage of using valence 
force field is that the potential energy distributions 
(PEDs) obtained here, in combination with the results 
of UB calculation’, can be used to identify the PEDs 
that are characteristic of the force field leading to the 
correct assignment of the vibrational modes as was 
done earlier®'® 

Method and results—The vibrational frequencies of 
DCA are taken from those reported by Devia and 
Carter’ me the revision suggested above. The method 
of Wilson’® is used to compute the force field. 

The molecule belongs to C, point group having 14 
vibrations of A’ species and 7 vibrations of A” species. 
In this note, A’ species are treated. The stretching and 
deformational modes of the methyl group belonging to 


Table 1—Final Force Constants of N,N-Dichloroacetamide 
(in units of mdyne/A, mdyne/rad and mdyne A/rad?) 


Symbol Coordinates Common Value 
involved atoms 
DIAGONAL CONSTANTS 
Stretch 
Kp C—H — 4.683 
Kp Cc—C — 3.850 
Kp C—N ro a ee 5.383 
K, N—Cl -- 2.782 
hs ; Cc=0 em 9.947 
Bend 
H, “HCH 
H, Lucc 
H, Fa 


N—Cl 0. 
int =P a. 
es 


_ 
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The assignment of the bands at 2936, 2857, 1427(A 
type asymmetric deformation of methyl group), 1001 
(C—C stretching mode) and 3001 cm * is independent 
of the force field used. That is, both the force fields agree 


A” type are also included. The structure parameters are 
the same as those employed by Devia and Carter’. The 
orthonormal set of symmetry coordinates used in the 
calculations can be had from the authors. 

The initial set of force constants, obtained by 
transferring from related systems, is refined using 17 
frequencies of DCA by the method of damped least- 
squares until convergence is obtained. The precautions 
necessary in this type of work as outlined in an earlier Obs. Calc. 


paper’ are also observed in this case. The final set of ie fren. 


Table 2—Potential Energy Distributions of 
N,N-Dichloroacetamide 


PED in internal force constants 


force constants thus obtained is given in Table 1. te a ; 
The potential energy distributions in internal force bee lot 
constants evaluated using well known expressions are 2936 2712 ih 
; : ; o 2857 2852 100Kp 
presented in Table 2. Contributions below 104, are not 1753 1744 71K,, 14H, 
shown. The vibrational assignments arrived at on the 1427 1435 85H, 
basis of PEDs in valence force field and PEDs in UB 1372 1388 62H,, 51H, 
force field are given in Table 3. The correct 1293 1293 57K p, 16K 
assignments, made on the basis of persistent PED on poe a ee 7K. 428 Rh 
. : - z e D> R> d> 2) 
contributions in the two force fields, are also given in 830 840 29K», 52K,, 15H, ‘ 
Table 3. eee . 582 601 18Kp, 17K,, 23H,, 10H, 
Discussion—The vibrational assignments presented 478 455 13K, 53Ky, 15H, 
in Table 3 are self-explanatory in that they give the 397 384 23K p, 45K4, 16H, 
assignment of vibrational frequencies of DCA in terms =e 5 14H, 115, 19K 
of mixing of vibrational modes. Yet a few explanatory = 2 ek 
A” species 


lines may be in order with regard to the field- 


dependence of vibrational frequencies. a ps ig 


1427-1432 82H, 
1036 —-:1026 90H, 


os eee 


Table 3— Vibrational Assignments of N,N-Dichloroacetamide 


Obs. Assignment on UB model Assignment on valence Si 
‘freq. (Devia & Carter’) force field Pecittaianc 
cm? (Present study) 
A’ Species 
aie es Vas (CH) 99 v,,(CH;) Vas(CH3) 
v,(CH;) 100 v,(CH;) vs (CH) 
i a i =O) + 12 v(C—N) 71 viC=O) WC a O) 
s(CH) + 12 »(CH,) 84 6,,(CH;) + 1 : . 
rik 82 6,(CH3) + 15 v(C—N) 91 Sete: st a gt ee 
64 ean, 32v(C—C) 53 v(C—N) + 18 v(C—C) (CN) + v(C—C) + 6(OCN 
+ 156(OCN) + 246,(CH,) + 110(OCN) 
are ss eee I5v(C—C) + 165(OCN) 60 »(CH) + 18v(C =O) »(CH3) 
; Ck )+ Llv(C—C) + 10d(OCN) 11 v(C—N) + 12v(C—C) v(C—N) + v(C—C) + 6(OCN) 
a0 ss INCE. 7d + 166(OCN)+ 13y(CH;) + y(CH;) 
pate + Se > SSC +3 ; Foon 2) 37 v,,(NCI,) + 28v(C—C) Vas(NCI,) + v(C—C) 
26v,(NCI,) (CCN+ 18 v(C—C)+ 336(CCN) + 16v,(NCI,) v(C—C) + 6(CCN) + v,(NCI,) 


478 Fo 5(OCN) + 22v,(NCI,) 
+ 37v,4,(NCI,) 
397-39 6(CCN) + 37¥,(NCI,) 


24 5(OCN) + 48¥44(NCI;) + 33v(C—C) 5(OCN) + v4,(NCI,) 


ben, 26 d(CCN) + 49v,(NCI,) + 23n(C—N) 65(CCN) + v,(NCI,) 
75 8(NCI,) + 12v,(NCI Oat 
20174 (NCI) + 165(OCN) + 145(CCN) 71 »(NCI,) + i1s0CN) HRCI + HOC 
2 
3001 100 v,,(CH,) scape : 
1427 «81-5 ,.(CH,, ‘Se 
es(CH) + 159(CH,) 91 5,,(CH het 
1036 85 y(CH,)+ 186,,(CH;) 91 (CH : : 
3 y(CHs) 


*1191 em! according to Devia and Carter! 


186 
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NOTES 


perfectly with regard to the assignment of the above 
bands. However, the two force fields give different 
assignments in the case of other vibrational modes as 
can be seen from Table 3. For example, the band at 478 
em”! is described as 6(OCN) + v s(NCI,) + v,,(NCI,) 
on the basis of UB calculations, whereas the valence 
force field describes it as 6(OCN) + Vas(NCI,) + v(C— 
C). In such cases, the assignment is based on persistent 
PED contributions in both the force fields. Hence the 
correct assignment for the above band is 5(OCN) 
+ Vas(NCI,). The contribution from v,(NCI,) in UB 
calculations and the contribution from v(C—C) in 
valence force field calculations are ignored treating 
them as characteristic of the particular force field. A 
similar procedure is adopted in assigning the other 
vibrational modes which can be found in the last 
column of Table 3. 


_ (i) The vibrations of the methyl group at 2936, 2857, 
1427, 1372, 3001, 1427 and 1036 are essentially pure. 
The rocking mode of this group at 1001 cm“! gets 
maximum PED contribution from )(CH3). 

Gi The amide I band at 1753 cm * is essentially due 

retching of C=O bond. The amide III band is 
Se sicassed by the C—N stretching character whereas 
estcanat of the amide IV band is distributed in the 
s near 1293, 1001 and 478 cm™*. 
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(ili) The bands around 397, 238 and 201 cm™! are 
primarily due to symmetric stretching, deformation 
and rocking modes of NCI,. The PED of asymmetric 
Stretching mode of NCI, appears in the bands at 830 
and 478 cm “!. 

The authors are grateful to Prof. K Venkata 
Ramiah, Member of the Union Public Service 
Commission, New Delhi, for helpful suggestions and to 
Prof. N A Narasimham, Head, Spectroscopy Division, 
Bhabha Atomic Research Centre, Bombay, for 
providing computational facilities. 
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CNDO/2 calculation is made for diformamide, C- 
fluorodiformamide and C-chlorodiformamide in frans-trans, trans- 
cis and cis-cis configurations. The charges, bond orders and dipole 
moments are discussed. 


Diformamide (DFA) is the simplest molecule 
containing the —CONHCO— group. Its study is of 
interest in relation to the structure of many 
biochemically important compounds in which the CO 
and NH bonds occur alternately. It has been 
established in the case of amides that no single valence 
bond structure could explain all their properties. This 
has been attributed to the delocalization of the 
carbonyl pi-electrons and lone pair electrons of 
nitrogen, resulting in the partial double bond character 
of the C—N bond?" °. In diamides, nitrogen lone pair 
and pi-electrons of two carbonyl bonds take part in 
conjugation. Halogen substitution results in additional 
pi-conjugation of lone pair electrons of halogen with 
pi-electrons of the amide system. Recently, we 
investigated the effect of such a conjugation in the case 
of trans-trans diacetamide (DAA) and _ trans-trans 
dipropionamide (DPA) by CNDO method®. It is the 
aim of this note to make similar Studies on 
diformamide (DFA), C-fluorodiformamide (FDFA) 
and C-Chlorodiformamide (CDFA) and to examine 
the effect of the above conjugation as well as the effect 
of halogenation at one of the functional carbon atoms 
in DFA on charge distributions, bond orders and 
dipole moments. 

Method of calculation—The CNDO/2 calcu- 
lations” of DFA, FDFA and CDFA in trans-trans, 
trans-cis and cis-cis confi guration are accomplished by 
programme QCPE 142 CNINDO®. The structure of 
the molecules is presented in Fig. 1. The structure 
parameters used are r(N—H)=1.0A, r(C—N) 
= 1402A, r(C=0)=1.240A, r(C—H)=1074 
r(O—F)=1.32A and r(C—Cl)=1.72A. The bond 
angles are assumed to be 120° each. The atomic 
orbitals used in the calculations are | for hydrogen; 2s 
2p for carbon, nitrogen, Oxygen and fluorine; and 3s, 3p 
and 3d for chlorine. \ 


Results and discussion—The net total charge (7 +9), 
net pi-charge and mobile bond orders of DFA, FDFA 
and CDFA are presented in Table 1. 

In DFAs the o-charge (this is total charge minus net 
pi charge) on hydrogen, fluorine and chlorine atoms at 
functional carbon are —0.027e.  —0.267e and 
—0.153e in trans-trans configuration as can be seen 
from Table 1. The net negative charge on the 
substituent means that the substituent withdraws 
electrons from the remaining molecule. Thus all 
substituents are o-withdrawing in nature. Their 
capacity to withdraw is in the order hydrogen 
< chlorine < fluorine for trans-trans configuration of 
DFA, FDFA and CDFA. The o-charges obtained for 
these molecules in trans-cis and cis-cis configurations 
are also, in general agreement with the above 
statement. 

The substituents behave differently regarding the 
conjugation through pi-bonding. The net charges in 
the pi-orbitals of substituents fluorine and chlorine are 
0.068e and 0.013e respectively in  trans-trans 
configuration. In trans-cis and cis-cis configurations 
this charge remains almost unchanged. Thus these two 
substitutents at the functional carbon are pi-donating 
type. Their capacity to donate is in the order chlorine 
< fluorine. 

Recently, Srinivas Rao® has studied some amides, 
thioamides and substitutions in them using CNDO/2 
method and concluded that the o-withdrawing effect of 
the substituents at functional carbon atom is in the 
order methyl group < amino group ~chlorine < 
fluorine; the pi-donating capacity of the substituents is 
in the order chlorine ~ methyl group < fluorine < 
amino group. Del Bene et al.'° arrived at similar 
conclusions from an ab initio study of some carbonyl 
compounds and substitutions in them using STO-3G 
basis set. The results obtained in the present work are 
in accordance with the conclusions of Srinivas Rao? 


(1) 


a1) 
Fig. |—Structure of: diformamide 


NOTES 


Net total charge (1 +0) 
N Co). 2° Of) 8 H(i) H N 
Diformamide —201 330 330 —266 -—266 -27 123 —24 252 
(TT) 
P Diformamide —205 330 330 —288 —288 -20 140 —20 254 
(TC) 
| Diformamide —208 325 325 —283 —283 jy Me By 17 256 
(CC) 
C-Fluoro- —222 334 532 —258 —313 -—199 140 14 241 
diformamide 
’ (TT) 
C-Fluoro- —229 330 540 —267 —333 —184 158 15 241 
diformamide | 
; (TO) 
; C-Fluoro- —231 330 529 -—275 —328 —195 174 4 244 
' (CC) 
C-Chloro- —193 331 388 -—260 —247 -140 137 15 251 
diformamide 
(TT) : 
C-Chloro- —200 332 385 —266 —270 -114 148 15 253 
, (TC) 
C-Chloro- —201 328 382 -—276 —263 —133 166 3 254 
diformamide 
(CC) 
' * Charges are in units of 10~* electronic charge. 


- *TT, TC and CC stand for trans-trans, trans-cis, and cis-cis configurations 
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and Del Bene et al.'!° Further, the same workers”’'® 
ad that o-effect predominates the pi-donating effect. 
n the present work also it is found that o-withdrawing 
Rae Seorine and chlorine when yingaiee 
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Table 1—Charges* & Mobile Bond Orders of the Molecules 


r whereas the C— 


Net n-charge Mobile bond order 


CC CAy.-6 OV) R N—C C=O N—C(1) C(1)=O(1) C(1)—R 
178 178 —305 —305 — 0.382 0.903 0.382 0.903 _ 

191 191 —318 —318 — 0.384 0.900 0.384 0.900 — 

185 185 —313 -—313 — 0.384 0.901 0.384 0.901 

173 207 —293 —396 68 0.372 0.907 0.378 0.852 0.286 
182 219 —302 —409 70 0.372 0.907 0.378 0.847 0.293 
183 214 —304 —406 68 0.373 0906 0.381 0.848 0.288 
173 182 —295 —324 13 0.377 0.907 0.378 0.864 0.217 
189 194 —314 —337 16 0.383 0.900 0.376 0.861 0.226 
183 190 —305 335 13 0.377 0.905 0.381 0.861 0.219 
respectively. 


configurations by 6.33 kcal/mol and 4.79 kcal/mol 


respectively. 
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Tanabe-Sugano diagrams are of great help to scientists working 
on the transition metal ions embedded in crystals in the optical 
absorption spectra. Since fixing the band positions with Racah 
parameter (C) involves a great amount of computer work, the cubic 
field energy matrices for d?, d?, d*, d°, d®°, d’, and d® electronic 
configurations are diagonalized. A few typical energy level diagrams 


mo _ drawn between Dq/B and E/B for certain values of C/B are 


— 


— 
ie 


- Presented. 


be In the faaatisic of the apical absorption spectra of 


‘transition metal ions embedded in crystals, we often 
use the well known Tanabe-Sugano diagrams!. These 
diagrams were drawn between E/B and Dq/B where E 


Is the nese? of the electronic level calculated by 


i the electrostatic and crystal-field 
‘2 B the electronic repulsion (Racah) 

and Dg the crystal-field splitting parameter. 
| agrams for d, 3, d*,d5,d°,d” and d®_ 


Two of the authors (YS) and (KP) are ies . 
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Lattice 


The generalized form of Born-Mayer equation without 
compressibility term has been modified and a polynomial of the type 


U, = Upuectk — 0.3729 x + 0.0695x? +...] 


has been suggested for rapid and easy calculation of lattice energy in. 


crystals. The relation is applied to calculate lattice eneegy in more 
than 40 different crystals. The results obtained show marked 
improvement over those reported by earlier workers. 


Recently, considerable efforts’ * have gone into the 
study of crystals and energy transformations in them 
for various processes. While studying crystals, one 
comes across the following generalized form of Born- 
Mayer equation for the calculation of lattice energy: 


U,= Use| -£| ) 
ro 


where 9 is a constant equal to 0.345 A, Ustec is the 
Madelung term given by: 


Km 
Usec =| ...(2) 


_. (Se ene 


Table 1—Values of Interionic Distances ro (A) and Lattice 
Energy U, (kJ mol~') for Alkali Halide Crystals 


Crystal rg UL Exptl** Born-Mayer 
(A) Present results results 
results 
LiCl 2.566 856.3 857.6 818 
NaCl 2.814 803.0 804.4 756 
KCl 3.139 717.4 719.1 687 
RbCl 3.285 689.9 691.7 661 
CsCl 3.560 673.6 675.5 621 
LiF 2.009 1003.0 1039.8 1000 
NaF 2.307 923.0 922.8 894 
KF 2.664 818.8 819.9 792 
RbF 2.814 783.4 784.8 756 
CsF 3.005 730.9 732.5 714 
LiBr 2.747 820.3 821.6 772 
NaBr 2.981 752.2 753.8 719 
KBr 3.293 688.8 690.6 os 
434 664.5 665.9 
ol : pages 598 
765.9 710 
705.3 670 
649.6 620 
628.7 600 
618.2 565 


Table 2—Values of Interionic Distance rp (A) and Lattice 
Energy U, (kJmol™') for Alkaline-Earth Chalcide and 
Heavy-Metal Halide Crystals 


Crystal rp U, Exptl** Born-Mayer 
(A) Present values results 
results 
MgSe 2.731 3305 3339 3103 
CaSe 2.962 3009 3038 2892 
SrSe 3.122 2918 2900 2764 
BaSe 3.302 2737 2761 2631 
MegTi 2.770 3116 3148 2879 
CaTe 3.179 2816 2841 2720 
SrTe 3.331 2767 2791 2611 
BaTe 3.500 2611 2632 2498 
BeO 1.649 4407 4540 4369 
MgO 2.105 3860 3890 3853 
SrO 2.580 3254 3332 3257 
NHF 2.630 740.5 738.1 752.1 
NHCl 3.34 684.3 681.2 656.7 
NHBr 3.51 656.4 649.4 630.2 
NHI 3.78 612.1 613.8 588.0 
TiF 2.59 848.2 828.4 811.6 
TiCl 3:33 686.0 698.7 658.4 
TiBr 3.44 666.4 689.2 639.7 
Til ee 04 636.6 673.8 611.3 
CuCl 2355 892.3 928.4 852.3 
CuBr 2.46 830.9 903.7 794.4 
AgF 2.46 910.9 886.6 847.6 
Cul 2.62 787.2 892.2 753.4 
AgCl pag | 800.3 849.4 766.5 
AgBr 2.89 771.1 824.3 738.9 
Agl 2.80 744.1 832.6 712.8 


*See Refs 8 and 9; **See Ref. 7. 


with k,,=NZ,Z2e?A where (Z,e), (Z2e) are the ionic 
charges, N the Avogadro number, A is the Madelung 
constant and rp is the interionic distance. From Eq. (1) 
it is evident that a straight line graph will be obtained 
between (Uetec = U,)/ Vevec and te: Thakur et al.4 
modified Eq. (1) to the following form: 


ad 


0 


U,= Ususe| 1085 — ..(3) 
Thakur et al.* also proposed a logarithmic type of 
lattice energy potential function. When the experimen- 
tal lattice energy values of different crystals are plotted 
against their corresponding ro values, it is found that 
this experimental curve does not depict a linear 
relationship between U, and ro but bears a curvature. 
The change in the values of U,; becomes very gradual 
for the crystals having larger values of ro. 

Further, the values of lattice energy for various 
crystals calculated from relations proposed by earlier 
workers do not agree with their corresponding 
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experimental values. In order to circumvent the above 
difficulties, we propose in this note the following 
empirical expression for the lattice energy function 


U, = Userec(k — ax + bx? — cx? + ...] ...(4) 


In this expression the constant k and the coefficients 
a, b, c etc. have been obtained iteratively using known 
standard data>~§ of U, and rg for a set of six to eight 
crystals. The final expression which is applicable to 
various crystals can thus be written as 


U, = Usrec[k — 0.3729 x + 0.0695 x? 
— 0.0086x? + 0.0008 x* + O(x°*)] — ...(5) 


where x = 1/rp and k = 1. The value of constant k in 
case of alkaline-earth chalcide crystals is 1.04. 

The above equation includes the combined effects of 
all weak mteractions (such as van der Waals, dipole- 
quadrupole term, etc.) existing between all unlike ion 
pairs. , 

From relation (5), lattice energy for more than 40 
different crystals, including diatomic and heavy-metal 
halide crystals, has been calculated. The results 
obtained are shown in Tables 1 and 2. 

For comparison we have also included in Tables 1 and 
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2 corresponding energy values obtained from Eq. (1) 
and those from experimental methods. It can be seen 
from these tables that the error in calculated lattice 
energy has been considerably reduced by using Eq. (5). 
The average percentage error in lattice energy 
predicted for alkali halide crystals is less than 0.18 %, 
while in the case of crystals shown in Table 2, the error 
is less than 0.8%. It is thus seen that the results 
obtained by us have a good agreement with the 
experimental values. 
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Tunnelling phenomenon for two different confining potentials 
depicting the leakage of scalar quarks has been discussed. The 
transmission coefficients for these cases have been evaluated 
numerically and compared with those for the linear potential. 


The concept of quarks was introduced by Gell-Mann’ 
and independently by Zweig? in 1964. Since then 
considerable efforts have gone in constructing 
potential models for quarks. Kang and Schnitzer* 
- constructed, from purely field theoretic considerations, 
a simple relativistic potential model to describe meson 
spectroscopy, using a linear function as the fourth 
component of a four-vector in Klein-Gordon equation. 
Kang and Schnitzer* have described mesons as bound- 
states of quarks and anti-quarks and their study has 
revealed the tunnelling and leakage of quarks. Iyer and 
Sharma*’* have studied quark confinement and the 
meson spectroscopy for a harmonic cum centrifugal 
potential in the relativistic limits. Ram®° and Ram and 
Halasa’ further discussed the tunnelling in terms of 
numerical results for linear and harmonic oscillator 
potentials respectively. They studied these potentials 
from purely phenomenological viewpoint. Further, it 
has also been observed that a potential providing 
confinement in the non-relativistic Schrodinger 
equation, gives rise to tunnelling’ ’ when used in 
relativistic equation. 
Motivated from the phenomenological viewpoint, 
we have studied, in this note, two confining potentials 
in the Klein-Gordon framework using WKB 
approximation. We have considered the tunnelling 
mer on by evaluating the values of the 
nission coefficients numerically. 
—The following potentials have been 


ve oh) 
- wath) 
h motion of a scalar quark moving in 
potential V;(r), its motion in the 


--(3) 


The radial form of Eq. (3) is 


d?X se ‘ 
apr + 2mLE - V§"(r)] X(r) =0 ...(4) 
X(r)=rRi(r), Wr) = Rir) YO, ¢) 

2 
f--—" 

2m 
and 

(l+1) E 2 

Vi"(r) = ( ) ~~ pag EP ee Fy 


ono mn” fy 2m? Se 


The chances of scalar quark leakage through the 
potential barrier depend on Vf"(r). When r = 0, V{"(r) 
= 2% and hence T =0. Again in Eq. (4) when !=0 
V(r) has a maximum value E2/2m at r= roe”’’. So at 
E =m, the potential barrier is only half as high as the 
rest mass m of the quark. Using WKB approximation 
and following Merzbacher®, one gets for transmission 
coefficient 


4 
f= 


1 2 
(20455) 
=exp| | Keodr | (7) 


Here r, and r, are the roots of the equation 


.. (6) 


K(r)=0 .. (8) 
and 
K (r) = [2mi Vi"(r) — E}]'? ...(9) 
Again considering potential V(r) one gets 

EB Eg 1 , 
Viti) =— gir +2 — 5 [air + eal” ...(10) 


For | =0, V#f(r) has a maximum value E7/2m at 


1/3 
- (* — $2 ) 
&1 

Results and discussion—We have obtained values of 
transmission coefficient 7 for lighter and heavier 
quarks moving with different energies through the 
potential barrier described by Eqs (5) and (10). For 
evaluating the integral occurring in Eq. (7), numerical 
method has been used. These values of T have been 
tabulated in Tables | and 2. For comparison 
corresponding values of T obtained for a linear 
potential® have also been included. 
Tt is seen from Table | .that from amongst the 
potentials considered in this work, the cubic 
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Table 1— Dependence of the Transmission Coefficient T for 
the 2 Potentials on the Quark Mass 


Quark j 
mass, — - — 
GeV Pot. (1) Pot. (2) Linear 
pot. 
(Ref. 6) 
0.1 0.61 0.63 0.63 
0.2 0.53 0.62 _ 
0.3 0.40 0.60 0.53 
0.5 0.14 0.55 O37 . 
1.0 46x 10~* 0.27 0.04 


The parameters chosen are: 2; = 1 GeV (Ref. 4); g. =0; E=1 GeV; 
ro = 1 (GeV) ' and g =1 GeV. 


interaction provides the maximum value for the 
transmission coefficient. It implies that chances of 
scalar quark leakage are more for this interaction when 
compared to other two potentials considered. Our 
observation clearly raises the question: Has this result 
been obtained because of applying Klein-Gordon 
equation or should this high value of T be interpreted 
as an indicator for the observability of light scalar 
quarks? The latter interpretation seems to be more 
reasonable for it leads to the existence of light scalar 
quarks as physical entities. 

Our above observation is consistent with that of 
earlier workers. For example Wilson’, in his gauge field 
theory has observed that the quarks are confined only 
in the strong coupling limit in which case the theory no 
longer remains Lorentz invariant. Similarly, 
Polyakov'® has also indicated that heating of 
confining gauge fields may lead to quark liberation. 
Further, La Rue et al.'' have claimed of having 
experimentally seen fractional charges in their 
superconducting levitation experiment involving 
niobium pellets. 
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Table 2—Variation of the Transmission Coefficient T with 
the Potential Parameters 


Quark T 
mass, ae 
GeV Pot. (1) Pot. (2) Linear 
ro = 1 ro = 0.95 Pot. 
g=1 g=2 g=!1 g,=0. g,=02 a=0.1 
0.10 0.62 0.63 0.62- 0.61 — 0.52 
0.15 0.61 0.62 0.61 = -- — 
0.20 0.59 0.61 0.58 0.55 0.59 = 
0.30 — — — 0.50 0.54 0.06 
0.50 0.28 0.48 0.27 0.36 0.43 4x10~* 
100 18x107* 0.10 — 0098 017 2x10" 
1.25 1.03x10~7’ 0.018 — 0,043 0.098 — 


1.35 642x10~'° 0.0076 — — —- as 


a ee Ee a ee a 


Table 2 depicts the variation in the value of 
transmission coefficient T with potential parameters. 
From Table 2 it is seen that the value of T for heavier 
quarks (m ~ 1 GeV) are relatively smaller. This may be 
ascribed to the reason that the heavier quarks move in 
the non-relativistic limit as (E—n)<m. This 
conjecture is in agreement with a similar observation 
made by earlier workers®. 
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